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INTRODUCTION 

Accurate  prediction  of  the  aerodynamic  flow  field  in  which  the 
helicopter  rotor  operates  is  critically  important  in  the  cal- 
culation of  helicopter  air  loads  and  blade  response.  In 
recent  times,  the  development  of  a model  for  nonuniform  inflow 
determination  was  carried  out  independently,  and  at  about  the 
same  time,  by  Piziali  and  DuWaldt  (ref.  1)  and  by  Miller  (ref. 
2).  The  many  applications  of  this  theory,  such  as  those  of 
references  3 and  4,  have  made  quite  clear  the  importance  of 
wake  structure  and  geometry  in  determining  blade  dynamics.  In 
particular,  the  outboard  portion  of  the  wake  of  each  blade, 
which  rapidly  rolls  up  into  a tip  vortex,  has  a significant 
influence  on  higher-harmonic  loading  (refs.  3 and  5)  and  over- 
all performance  (ref.  6). 

Since  the  wake  has  such  a strong  effect  on  the  blade  loading, 
it  will  obviously  strongly  influence  the  performance  char- 
acteristics of  a rotor.  Tnis  is  especially  true  if  (1)  two 
overlapping  rotors  exist  in  the  system  so  that  the  wake  from 
one  will  pass  through  the  other  rotor,  (2)  the  flight  con- 
dition is  such  that  the  wake  remains  nearby,  or  (3)  large 
radial  and/or  azimuthal  variations  occur  in  the  loading  dis- 
tribution so  that  strong  vortices  are  shed  into  the  wake.  The 
effect  of  the  wake  on  the  stability  and  control  derivatives  is 
not  quite  so  obvious.  A system  which  would  appear  tc  be 
stable  if  uniform  inflow  is  assumed  in  analyzing  the  loads 
might  turn  out  to  be  less  stable  (or  even  unstable)  if  the 
actual  wake-induced  velocity  distribution  were  used  in  the 
analysis. 

Because  of  the  importance  of  considering  an  accurate  wake  re- 
presentation, especially  for  advanced  configurations,  a study 
has  been  conducted  here  tc  investigate  the  effects  of  the 
wake  on  rotor  performance  and  stability.  The  basic  wake  geom- 
etry and  blade  loads  programs  used  here  are  described  in 
reference  7. 


Methods  for  predicting  the  rotor's  aerodynamic  flow  field  have 
recently  undergone  extensive  refinement,  primarily  due  to  the 
need  for  improved  predictions  and  the  availability  of  high- 
speed, large-capacity  digital  computers.  These  refinements 
have  resulted  in  more  realistically  detailed  models  for  calcu- 
lation of  nonuniform  wake-induced  inflow.  Some  of  these  models 
are  discussed  in  references  7 through  11.  Corresponding 
blade  air  load  and  blade  response  models  have  been  developed, 
so  compatible  blade  loading  and  response  models  can  be  used. 

The  improved  models  for  the  wake-flow  prediction  have  provided 
a good  qualitative  approximation  to  average  measured  wake  flows 
(refs.  7,  12). 


The  determination  of  wake- induced  flow  at  the  rotor  blades  has 
been  the  most  difficult  part  cx  . ,e  model  to  develop  satis- 
factorily. This  is  pvir.f.ril**  a result  of  the  sensitivity  of 
wake-induced  flow  to  the  magn ix ude  of  the  relatively  small  dis- 
tances between  the  tip  vortex,  of  one  blade  and  the  subsequent 
blade,  esp.eci  ily  on  the  advancing  side.  Similarly,  because 
of  the  interactions  o‘  -uo  various  wake  elements,  t.ie  location 
of  the  lii’.-c  -.rd  wake  s*  n-M-ure  also  has  a significant  effect 
ot;  blade  h.  Tnus  a reasonably  accurate  definition  of 

the  waki  geometry  is  needed  in  ord.:r  that  wake  effects  can  be 
proper Ly  evaluated-, 

A rigor-)i  ? cV.:ermiria>  • vi;  or  *.ne  structure  and  motion  of  the 
wake  of  r.  .rocor  is  a4*  ex  ■ * ly  arduous  task.  The  flow  is  un- 
steady and  t'lreu- nimer.sion-tl,  and  allows  no  convenient  basis 
for  linearisation.  As  a result,  in  the  past  the  wake  has  been 

assumed  to  be  a fixed  skewed 'helical  pattern  for  forward 
flight  conditions,  and  wake  distortions  have  generally  been 
neglected.  Studies  reported  in  refs.  13  and  14  indicated  the 
possibility  of  predicting  wake  distortion  effects  by  consider- 
ing discrete  elements  of  the  tip  vortex  in  the  wake  of  a 
propeller  or  helicopter  rotor. 

In  reference  13  the  flow  due  to  a rotor  hovering  near  a ground 
plane  was  approximated  by  axinymmetric  flow  due  to  a succession 
of  vertex  rings  released  from  the  rotor  plane.  A fairly  good 
approximation  to  the  physical  flow  being  modeled  was  obtained. 
Near  the  rotor  plane  the  interior  flow  was  poorly  predicted; 
it  is  believed  that  this  resulted  because  the  inboard  wake 
structure  was  not  represented.  In  reference  14,  the  flow  due 
to  a rotor  in  forward  flight  was  calculated.  The  wake  of  each 
blade  was  represented  by  a single  concentrated  tip  vortex 
having  a finite  core  of  rotational  fluid.  Good  agreement  with 
flow  measurements  some  distance  from  the  rotor  plane  was  ob- 
tained. Again,  however,  it  is  believed  that  omission  of  the 
inboard  wake  structure  is  the  reason  that  large  errors  in  the 
predicted  flow  in  or  near  the  rotor  plane  were  noted. 

In  order  to  check  out  the  wake  geometry  analysis  developed  by 
Sadler  in  reference  7,  comparisons  were  made  with  the  induced 
velocity  distributions  in  the  wake  of  a model  rotor  as  mea- 
sured by  Heyson  and  Katzoff  in  reference  12.  The  self-deform- 
ing wake  which  was  developed  by  Crimi  in  reference  14  and 
involved  only  the  tip  vortex  model  of  reference  14  produced 
good  results  outside  of  and  at  the  edge  of  the  wake  but  pro- 
vided poor  comparisons  with  the  measured  values  over  the  inner 
portion  of  the  wake,  especially  up  near  the  rotor.  The  full- 
mesh  model  of  reference  7 generally  compared  well  with  the 
measured  values  throughout  the  entire  region.  These  improved 
induced-velocity  calculations  along  the  whole  blade  will 
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provide  a better  basis  for  computing  the  loading  on  the  rotor. 
This  full-mesh  model  then  is  the  one  used  when  determining 
wake  effects  for  the  cases  conside*'2d  in  this  report.  Some 
attempts  have  been  made  to  account  for  wake  distortion  in 
blade  load  analysis.  In  the  calculation  of  VTOL  propeller  per- 
formance in  reference  15,  determination  of  wake  distortion  was 
made  a part  of  the  overall  calculation.  In  this  wake,  consid- 
eration was  limited  to  the  hover  case  which  made  the  flow 
steady  relative  to  the  blade.  In  a similar  approach,  the  the- 
oretical results  of  reference  13  were  used  in  reference  6 to 
, make  an  estimate  of  wake  contraction  for  VTOL  propellers  in 
yawed  flight.  Although  the  procedures  used  in  these  studies 
are  limited  to  rather  special  cases,  they  do  indicate  the  im- 
portance of  accounting  for  wake  distortions  in  the  performance 
analysis  of  rotors  and  propellers. 

The  present  program  has  the  capability  of  predicting  wake  geom- 
etry and  wake  flow,  including  the  induced  flow  in  the  rotor 
plane  for  use  in  blade  loads  calculations,  for  a system  of  one 
or  two  helicopter  rotors.  If  two  rotors  are  used  for  the  heli- 
copter model,  they  may  be  arbitrarily  located  with  respect  to 
one  another,  have  blades  of  different  properties,  and  rotate 
in  the  same  or  different  directions;  and  the  wakes  from  all 
blades  are  allowed  to  interact  and  to  affect  the  induced 
velocity  distribution  on  each  of  the  rotor  blades.  The  blade 
loads  program  uses  input  from  both  an  independent  blade  fre- 
quency and  mode  shape  prediction  program  and  from  the  wake 
geometry  program  discussed  herein  to  calculate  blade  loads  and 
blade  response.  The  forces  and  moments  transferred  to  the  hub 
by  all  of  the  blades  are  then  computed  in  order  to  find  the 
; values  for  the  performance  parameters.  Rerunning  the  blade 

loads  and  response  program  with  values  for  the  blade  pitch 
I settings  or  aircraft  attitude  which  are  perturbed  from  the 

; trim  settings  then  produces  the  stability  and  control  deriva- 

tives. 

: 

s The  computer  programs  which  were  used  for  this  investigation 

t were  modifications  of  those  developed  in  reference  7.  The 

l major  modifications  involved  were  (1)  the  calculation  of  the 

\ forces  and  moments  transmitted  to  the  hub,  (2)  inclusion  of 

; program  controls  for  computing  static  stability  and  control 

| derivatives,  (3)  use  of  radial  variations  for  the  airfoil  type, 

| chordlength,  and  nonlinear  twist  as  required  by  the  advancing 

| ' blade  concept  rotor  (see  ref.  16),  and  (4)  inclusion  of  air- 

foil flap  deflection  effects  in  the  calculation  of  aerodynamic 
'<  loads  for  the  wake  geometry  and  blade  loads  programs  and 

t elastic  twist  for  the  wake  geometry  as  was  required  for  the 

• controllable  twist  rotor  (see  ref.  17). 
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ANALYSIS 


A rotor  performance  analysis  requires  that  the  forces  and 
moments  which  are  transferred  to  the  hub  be  calculated.  These 
forces  and  moments  are  the  result  of  the  loading  acting  on  the 
individual  blades  of  the  rotor.  If  the  blades  are  flexible, 
then  the  loading  will  depend  on  the  response  of  the  blades  to 
that  loading,  and  so  iterations  must  be  carried  out  on  the 
blade  loads  and  responses  until  they  become  compatible  with 
each  other.  The  aerodynamic  loading  will  be  dependent  upon, 
among  other  things,  the  velocity  distribution  which  is  induced 
by  the  wake  from  the  rotor.  Thus,  a full  analysis  of  the  per- 
formance parameters  and  stability  and  control  derivatives 
would  require  that  (1)  the  wake  geometry  and  wake  influence 
coefficients  be  determined  so  that  the  wake-induced  velocities 
can  be  calculated,  (2)  an  interactive  blade  loads  and  response 
computation  be  carried  out,  and  (3)  upon  convergence  of 
the  blade  loads  and  response  calculations,  the  forces  and 
moments  which  are  transferred  to  the  hub  be  calculated 
using  these  final  values  for  the  loading  distribution  on  the 
individual  blades.  Perturbing  the  aircraft  attitude  or  blade 
pitch  settings  and  rerunning  the  blade  loads  and  response  cal- 
culations woull  then  provide  new  values  for  the  performance 
parameters,  and  these  new  values  would  be  used  to  determine 
the  stability  and  control  derivatives. 

The  overall  model  and  program  arrangement  is  shown  in  Figure  1. 
The  model  and  problem  formulation  may  be  conveniently  thought 
of  in  three  sets:  one  dealing  with  the  wake  model  and  associated 
calculations  , the  second  dealing  with  determination  of  the 
aerodynamic  blade  loads  and  response,  and  the  last  section 
concerned  with  computing  rotor  performance  and  stability  and 
control  derivatives.  The  combination  of  these  analyses  and 
corresponding  computer  programs  arc  directed  toward  the  pre- 
diction of  blade  loads  and  stability  and  control  derivatives 
for  helicopters  in  steady  forward  flight  or  hover  in  which  the 
effects  of  a free  wake  and  flexible  blades  are  included.  One 
or  two  rotors  may  be  modeled,  and  if  two  rotors  are  used,  the 
blades  on  one  rotor  may  differ  physically  from  the  blades  on 
the  other  rotor.  Blade-wake  interactions  are  allowed,  but  no 
mass  or  elastic  couplings  between  blades  on  a rotor  or  between 
blades  on  different  rotors  are  allowed.  While  the  number  of 
blades  per  rotor  is  arbitrary  {subject  to  practical  limits  of 
the  dimensioned  variables  in  the  computer  program) , each  rotor 
is  assumed  to  have  the  same  number  of  blades  and  to  have  the 
same  rotor  speed. 


4 


,k-  -- 


..  





WAKE  GEOMETRY  MODEL  AND  FORMULATION 


The  wake  geometry  is  calculated  by  carrying  out  a process  simi- 
lar to  the  startup  of  a rotor  in  a free  stream.  The  blades 
are  located  at  specified  azimuthal  and  flapping  positions, 
without  any  wake  vortices.  The  blades  then  rotate  through  an 
azimuthal  increment,  At|»,  and  shed  and  trail  vortex  elements  of 
unknown  strength,  but  with  known  positions.  The  strengths  of 
the  vortices  that  are  shed  immediately  behind  the  blade  are 
then  determined,  and  include  the  effect  of  their  own  self- 
induced  velocities.  An  estimate  of  the  blade  loads  that 
result  is  then  determined,  without  the  effects  of  blade  flexi- 
bility being  included.  All  vortex  element  end  points  not 
attached  to  the  blade  are  then  allowed  to  translate  as  the 
blade  is  stepped  forward  for  a time  At  = At p/n,  where  n is  the 
rotational  speed.  This  completes  a typical  first  step  in  the 
wake  geometry  calculation. 


Subsequent  steps  are  similar.  In  this  manner  arrays  of  dis- 
crete shed  and  trailing  vortices  are  generated  immediately 
behind  the  blades  with  strengths  which  correspond  to  approxi- 
mate blade  loads.  These  arrays  have  stepwise  radial  and 
azimuthal  strength  variations  so  that  total  circulation  is 
conserved.  The  arrays  of  shed  and  trailing  vortices  which  are 
generated  immediately  behind  the  blades  are  referred  to  as  the 
full  mesh  wake.  Comparisons  of  wake  flows  predicted  using 
this  wake  model  for  the  entire  wake  with  experimental  measure- 
ments indicated  that  retention  of  shed  elements  with  a coarse 
mesh  resulted  in  poor  induced  velocity  predictions,  and  that 
use  o!  a fine  mesh  increased  running  time  to  an  unacceptable 
level.  Therefore,  the  full  mesh  wake  was  used  to  represent  the 
wake  immediately  behind  the  blades,  and  a modified  wake  model 
was  developed  and  implemented  for  use  in  the  representation  of 
the  remainder  of  the  wake,  as  shown  in  Figure  2. 


The  modified  wake  consists  of  trailing  vortices  only,  so  vor- 
ticity  is  not  conserved.  The  wake- induced  velocities,  wake 
distortions,  and  other  calculations  are  essentially  the  same 
for  both  the  full  mesh  wake  and  for  the  modified  wake  portions 
of  the  wake  model. 


Determination  of  Wake  Geometry 


The  right-handed  Cartesian  coordinate  system  used  in  the  for- 
mulation for  the  wake  geometry  calculations  has  the  x-axis  in 
the  downstream  direction  {positive  toward  the  rear  of  the 
helicopter,  ^=0)  and  the  z-axis  vertical  and  positive  up,  as 
indicated  in  Figure  3.  A second  rotor,  if  one  is  used  in  the 
model,  may  be  located  arbitrarily  within  this  system.  Given 
the  flight  conditions  and  appropriate  rotor  and  control 
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parameters,  wake  geometry  calculations  are  essentially  those 
required  to  compute  vortex  element  strengths,  vortex  induced 
flowr  and  motion  of  the  vortex  elements. 

The  Biot-Sava:^  i*v  is  the  basic  relationship  used  for  calcu- 
lating vertex  > -meed  flows,  and  gives  the  fluid  velocity  g.  at 
a point  located  !./  the  vector  s due  to  a vortex  with 

circulation  r 


as 


n | 


r (s ) s xd£ 
77  I 


1 


where 


£l 


= s - s 
-P 


In  order  to  compute  total  vortex  induced  flow  at  a point,  the 
integral  is  taken  over  all  vortex  elements  in  the  flow.  For  a 
straight  element,  as  shown  in  Figure  4,  with  end  points  lo- 
cated at  A and  B,  the  induced  velocity,  q , at  C is  given  by 


% a 773  <cos  eA  - cos9b) 


where  r is  the  strength  of  the  element  between  A and  B.  When 

computing  the  vortex  induced  flow  at  a point  due  to  an  ad- 
joining vortex  element,  the  preceding  relation  becomes 
indeterminate.  The  calculation  of  the  induced  flow  in  this 
case  is  discussed  in  reference  7. 

The  blade  circulations  are  calculated  as  follows.  The  veloc- 
ities V and  U,  normal  and  tangential  to  the  plane  normal  to  the 
shaft  axis,  are  given  approximately  by 

v * VB.  “ " w (l) 

U = nr  + sin  <J>  (2) 

where  a is  the  shaft  axis  angle  with  respect  to  the  free 
s 

stream, positive  aft,  a&  is  the  forward  tilt  of  the  rotor  plane 

with  respect  to  the  shaft  axis  due  to  flapping,  and  w is  the 

induced  downwash  due  to  the  wake  (w=0  at  startup,  i.e., 
there  is  no  downwash  at  the  startup) . 


i 
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Using  small  angle  approximations,  the  lift  per  unit  span  is 
given  by 

l = l/2pU2cc^  (3) 

where  the  chord  length,  c,  may  change  with  radius.  For  linear 
aerodynamics,  the  lift  coefficient,  c^,  can  be  expressed  as 

Ct“  c*  (a  - “L  > (4> 

P.  0 

where  c^  is  the  lift  curve  slope  and  may  change  with  radius 
a 

due  to  a radial  airfoil  distribution,  o is  the  aerodynamic 

angle  of  attack,  and  aT  is  the  angle  of  attack  for  zero  lift. 

Jb 

0 

The  aerodynamic  angle  of  attack  is  the  stun  of  the  geometric 
pitch,  otg,  and  the  angle  due  to  the  normal  and  tangential 

velocities  so  that 

a * a + tan"1  (V/U)  (5) 

y 

The  geometric  pitch  is  the  sum  of  the  collective  and  cyclic 
angles,  the  built-in  twist,  $ (r) , and  the  elastic  twist, 

♦ (r,i|>),  so  that 

o*0  - A.  costp  - B.siniJ/  + <|>(r)  + 4> 

y 0 X X 6 

The  elastic  twist  can  be  expressed  as 

N 

* a (r)  +JL  [a  (r)  cos  n\J>  + b (r)  sin  n<J»l 
6 o n— i [ n n j 

where  a^  and  b are  the  harmonic  coefficients  of  the  elastic 
n n 

twist  with  respect  to  if>. 

The  zero  lift  angle  is  essentially  the  result  of  camber  in  the 
airfoil.  A flap  deflection  produces  an  effective  camber,  and 
thus  the  zero  lift  angle  will  change  with  flap  deflection.  A 
good  representation  of  the  zero  lift  angle  is 


a.  (0)  + 
Lo 
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The  terms  a^  (0)  and  da^  /d6^  are  illustrated  in  Figure  5. 
o o 1‘ 

Both  aL  (0)  and  da^  /d6^  are  usually  negative  or  zero, 
o o 

The  Kutta-Joukowski  law  states  that 

i = purfa  (6) 

where  is  the  bound  circulation.  Combining  the  various  re- 
lations given  in  Equations  (1)  through  (6)  and  assuming  that 
V/U  is  very  small  yields  a bound  circulation  of 


= l/2cc 

a 


(a  -o_  ) (nr  + V-.sin'l') 

g L i 

v o 


+ Vs  * V 


-w 


(7) 

For  a set  of  rotor  blades  which  have  stepwise  radial  and  azi- 
muthal circulation  variations,  the  above  equations  may  be 
thought  of  as  applying  to  each  radial  and  azimuthal  location 
independently.  The  wake- induced  velocity  on  the  blade,  w,  is 
made  up  of  velocities  due  to  known  circulations  in  the  wake 
and  to  unknown  circulations  at  the  blade,  and  may  be  written 
in  the  form 

w<ri'*k>  = Wk>  + l | “ij'WVV  (81 


where  w^r^,^)  is  the  induced  velocity  due  to  all  known  wake 
circulations,  r (r 


is  the  blade  circulation  at  r£ 


, y and 


a ^ j <r i * ^k^  an  influence  coefficient  which,  when  multiplied 

by  the  circulation  r (r  gives  the  induced  velocity  of  that 

JC  J 


element  at  * T^e  sumrnat^ons  °ver  the  indices  i and  j in- 

dicate a summation  over  all  radial  sections  of  all  blades  at 
their  respective  azimuthal  positions.  Then  a set  of  equations 
for  all  r's  may  be  obtained  from  Equations  (7)  and  (8)  and  is 
of  the  form 


ik  = 3 cla 


'-l  l °ijik 

. j 


(y  aL  )(nri  + Vf*sin  \)  + Vf(ag  - ap>  - w^ 


(9) 
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Here  r^  is  equivalent  to  r(r^,  tl^)  and  occurs  on  both  sides 

of  the  equation.  This  equation  is  solved  with  a simple  iter- 
ative procedure.  The  procedure  is  as  follows.  The  terms  on 
the  right-hand  side  of  Equation  (9)  which  do  not  involve  the 
r’s  are  used  as  an  initial  estimate  for  the  r^'s.  Then  the 

r^^'s  are  evaluated  from  the  equation,  and  as  each  is  evalu- 
ated, it  replaces  the  appropriate  r,.  {on  the  right-hand  side) 

ij 

Once  a complete  set  of  r^'s  is  determined,  it  is  compared  to 

the  previous  set,  and  convergence  is  assumed  when  the  sum  of 
the  squares  of  the  differences  of  the  sets  of  r's  is  less 
than  0.00005. 


After  the  blade  circulation  values  are  determined  for  a 
given  azimuthal  location  of  the  blades,  vortex- induced 
velocities  are  computed  at  all  vortex  element  end  points  in 
the  wake.  The  blades  are  then  advanced  through  an  azimuthal 
increment,  cup,  and  each  vortex  end  point  which  is  not 
attached  to  the  blades  is  allowed  to  translate  for  the  time 
period  Aiji/ft,  with  the  resultant  velocity  due  to  the  free 
stream  and  induced  velocities.  The  entire  computational  pro- 
cess for  the  new  rotor  blade  positions  is  then  repeated. 

That  is,  new  wake-induced  and  free-stream  velocities  at  the 
blades  are  computed  and  used  to  determine  new  blade  circu- 
lations. Then  new  vortex  end  point  velocities  are  computed, 
etc.  Each  blade  advance  results  in  an  additional  set  of  shed 
and  trailed  vortices  being  added  to  the  wake.  The  number  of 
revolutions  of  wake  retained  for  actual  computational  pur- 
poses is  restricted  by  an  input  to  the  computer  program. 

Wake  Flow  and  Wake-Induced  Velocity  Influence  Coefficient 
Calculations 


Wake  flow  is  calculated  by  using  the  same  basic  programming 
which  computes  flow  at  a point  due  to  an  arbitrarily  located 
vertex  element,  excepu  that  the  position  of  the  point  is 
specified  by  input,  and  the  flow  is  averaged  and  nondimen- 
sionalized.  Input  for  the  blade  loads  part  of  program  is 
calculated  by  manipulating  and  properly  subscripting  numbers 
equivalent  to  the  a ^ j (r^,^)  and  similar  numbers  used  in  the 

computation  of  wR.  These  numbers  are  used  in  the  solution 

of  an  equation  in  the  blade  loads  program  which  is  approxi- 
mately equivalent  to  Equation  (9).  Both  the  wake  flow 
calculations  and  the  evaluation  of  the  input  to  the  blade 
loads  program  are  computed  only  after  a specified  number  of 
revolutions  have  been  done.  Thus,  the  input  (including  both 
o and  r type  quantities  as  discussed  above)  from  the  wake 
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geometry  program  to  the  blade  loads  program  is  based  on 
approximate  specified  blade  motions  and  an  approximately  re- 
petitive wake.  The  repetitive  nature  of  wake-induced  effects 
has  been  determined  by  visual  checking  of  data  to  occur 
after  approximately  nR/V^ir  rotor  revolutions. 

BLADE  LOADS  AND  RESPONSE  MODEL  FORMULATION 

The  right-handed  coordinate  system  used  in  the  calculation  of 
(periodic)  blade  circulations  and  blade  response  is  located 
such  that  the  z-axis  is  fixed  to  the  shaft,  directly  upward, 
the  x-axis  is  downstream,  blade  azimuth  angle,  \p , is  mea- 
sured with  respect  to  the  x-axis,  y is  positive  on  the 
advancing  side,  and  the  distance  radially  outward  from  the 
axis  of  rotation  on  a given  blade  is  denoted  by  r. 

Blade  Natural  Frequency  and  Mode  Shape  Determination 


The  program  which  calculates  blade  natural  frequencies  and 
mode  shapes  is  an  independent  program  which  was  developed  by 
RASA  and  could  be  replaced  by  any  approximately  equivalent 
program.  The  necessary  output  for  any  such  program  in  order 
that  it  might  be  used  as  input  for  the  blade  loads  program 
is  the  natural  frequencies  together  with  the  corresponding 
mode  shape  quantities,  i.e.,  flapping  and  edgewise  displace- 
ments, slopes,  shear  forces  and  moments,  and  torsional 
deflection  and  torque.  These  mode  shape  quantities  are  to  be 
defined  at  the  location  of  the  point  masses  of  the  lumped 
parameter  model,  with  mode  shape  magnitudes  adjusted  to  give 
unity  generalized  mass.  The  lumped  parameter  model  lengths, 
masses  and  inertias,  mass  eccentricities,  offsets  of  elastic 
axis  from  pitch  axis  and  midchord,  and  twist  distribution 
(but  not  modal  bending  or  torsional  stiffnesses)  are  used  in 
the  blade-response  program. 

The  blade  frequency  program  developed  by  RASA  is  described 
in  references  18  and  19,  and  the  model  used,  development  of 
input  data,  and  operation  of  the  blade  frequency  program  are 
discussed  in  more  detail  in  those  reports.  Briefly,  the 
model  used  for  the  real  blade  is  a lumped  parameter  approxi- 
mation consisting  of  uniform  massless  elastic  beam  sections 
under  tension  due  to  centrifugal  loads,  with  point  masses  and 
inertias  located  at  the  ends  of  the  massless  lengths.  A 
modified  transfer  matrix  approach  is  used  in  determining  the 
natural  frequencies  and  mode  shapes.  (See,  for  example, 
ref.  20.)  The  natural  frequencies  and  corresponding  mode 
shape  quantities  are  used  in  the  calculation  of  the  response 
of  the  flexible  blades  to  aerodynamic  and  inertia  loads  in 
the  blade  loads  and  response  part  of  the  program. 
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Determination  of  Aerodynamic  Loading 

The  aerodynamic  loading  at  a given  radial  and  azimuthal  sta- 
tion is  derived  from  the  flow  experienced  by  the  blade 
section,  as  sketched  in  Figure  6.  The  geometric  incidence  of 
the  section  with  respect  to  the  rotor  plane  is  the  sum  of  the 
rigid-body  pitch  angle  0(i}»),  blade-  twist  <Hr) , and  torsional 
deflection  $e(r,i|i).  The  velocity  component  tangent  to  the 

rotor  plane,  U, is  given  by 

U * fir  + V,  cos  o simp  (10) 

f s 


i 

» 

j 
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and  that  normal  to  the  rotor  plane,  V, is  given  by 


V * V-  sin  a + In (r , 4»)  - w - w - V-  cosa  cosijisinB 

i S CIS 

(ID 

vrtiere  A is  the  plunging  velocity  of  the  section  due  to  the 
response  of  the  blade,  and  w is  the  chordwise  average  wake- 
induced  downwash,  and  is  given  by 


wtri'*j>  ‘ l l ‘W’VV  r(V*n) 

J m n j 


where  a (r.  ,i|;.)  is  the  wake-induced  downwash  at  (r.,<K)  due 

mn  ID  ID 

to  unity  bound  circulation  for  blades  which  were  located  at 

(rm»^n^  ^rm'^n+N  etc*'  we  is  the  climb  rate,  and  B is  the 
S 

blade  flapping  angle  relative  to  the  shaft. 


As  is  indicated  in  Figure  6,  the  airfoil  section  is  re- 
placed by  a vortex  distribution  of  strength  y(r,5,t|i)  along 
the  chord.  This  distribution  is  adjusted  to  make  the  flow 
at  the  section  tangent  to  the  chord,  which  relates  r to  the 
w and  yields  the  basic  relationships  governing  the  aerody- 
namic loading.  More  details  of  the  derivation  of  r are 
given  in  reference  7.  The  resulting  equation  for  r is 
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blade  semichord  at  r^,  feet 

lift  coefficient 

2 2 1/2 
(U  4V)' 

9<*J+1)  - «<*j.1)  + ♦e<ri'*j+1)  - ♦e(ri'*j-l) 

sum  of  collective  and  cyclic  pitch  angles, 
radians 

torsional  blade  elastic  deflection,  radians 


In  equation  (12a)  the  lift  coefficient,  c^,  is  a function  of 

angle  of  attack,  o,  and  Mach  number,  M,  as  defined  by  aero- 
dynamic coefficient  subroutines.  The  total  angle  of  attack, 
a,  is  defined  by 

a (ri#ipj)  = eu.)  4 ♦eCri,*^)  4-  ^(ri)  4 tan  *Vi;./Ui;.) 

(12b) 

Equations  (12a)  and  (12b)  represent  a set  of  nonlinear 
equations  in  the  strengths,  r.  The  nonlinearity  is  a result 
of  nonlinear  aerodynamic  coefficient  definitions,  and  of  the 
nonlinear  dependence  of  r upon  itself  (as  contained  in  w and 
its  contributions  to  V,  then  o,  and  finally  c^) . Solution 

of  Equations  (12a)  and  (12b)  is  therefore  accomplished  in  an 
iterative  manner,  and  is  discussed  in  detail  in  reference  7. 


The  lift,  drag,  and  moment  per  unit  span  are  readily  calcu- 
lated once  blade  circulations  have  been  obtained.  Resolving 
the  forces  into  components  normal  and  tangential  to  the  rotor 
plane,  F and  F , respectively,  results  in  the  following  ex- 

pressions: 


F = pbA 
zij  1 

0bA 

+ JZt 


U(UC|  4 Vcd) 


j+X  - j+l  + 2bi  (omu,1+l  ' (cmc,i-l 


Fx£j  ’ “bi[0(-Vc1  + Mij 
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(13) 

(14) 


5 * 7iiT««irTT 


' T^ r*^Aw  ^'*sas7*c 


where  c,  » section  drag  coefficient  and  c = section  moment 
a m 

coefficient  (about  midchord).  The  moment  about  midchord,  Mq, 
is  given  by 


M0(ri,^;j)  « pb. 


2(ciru  }ij  " ^E(rj+1  * r^“l)i 


3irb. 


8 (At) 


a - 2a  + a 

9*j.i  *?.»  a -i 


*j+l 


‘j  'j“l 


+ (5^i  ~ S,_,)/2 


j+1  ”j-l' 


(15) 


where  K - local  blade  spanwise  slope,  radians. 


Calculation  of  Aerodynamic  Coefficients 

Three  types  of  procedures  are  available  in  the  blade  loads - 
response  program  for  determining  the  lift,  drag  and  pitching 
moment  coefficients  at  the  midchord  of  the  airfoil.  The 
first  type  uses  polynomial  curve  fit  techniques,  the  second 
type  interpolates  linearly  on  angle  of  attack,  and  the  third 
type  performs  a triple-interpolation  for  angle  of  attack, 

Mach  number,  and  flap  deflection. 

Curve  fit  techniques  have  been  applied  to  the  lift,  drag,  and 
pitching  moment  data  for  the  NACA  0012  and  NACA  0015  in  order 
to  determine  the  coefficients  of  polynomials  which  would 
closely  fit  this  data  over  various  ranges  of  the  angle  of 
attack.  The  data  (and  consequently  the  polynomials)  is  for 
incompressible  aerodynamics  and  covers  all  angles  of  attack 
from  -180°  to  +180°.  Corrections  for  compressibility  are 

made  by  dividing  all  terms  by  J 1-M2 . Further  details  are 
given  in  references  7 and  21  for  the  NACA  0012  airfoil,  and 
in  references  22  and  23  for  the  NACA  0015  airfoil. 

The  airfoil  shapes  used  on  the  model  for  the  wind  tunnel 
tests  for  the  advancing  blade  concept  rotor  ranged  from  an 
NACA  0006  at  the  tip  to  an  NACA  0030  at  the  root.  Among  the 


13 


airfoils  tested  and  reported  in  reference  24  are  the  NACA 
0006,  0009,  0012,  0015,  0018,  0021,  and  0025  airfo\ls.  The 
curves  for  the  sectional  lift,  drag  and  pitching  moment  co- 
efficients that  are  given  in  reference  24  were  ured  to  compile 
tables  which  could  be  used  for  the  ABC  rotor.  Using  a linear 
interpolation  to  find  the  lift  coefficient,  for  example,  at 
o,  which  lies  between  and  cu  we  find  that 


c*(o)  ' c£(ai)  + da 


cl(ai+l)  " ct{ 


ai) 


= Cl-da)  ct(a.)  + daCjl(ai+1) 


(16) 


where 


a ~ ai 

ai+rai 


The  data  is  for  incompressible  tests,  and  compressibility  was 
included  through  the  relation 


c, (a,M) 

A» 


Ct(a)/ 


Values  for  c.  and  c were  obtained  in  the  same  manner  as  des- 
a m 

cribed  above. 

The  tests  in  reference  24  were  conducted  at  angles  of  attack 
from  -2°  to  28o.  The  angles  of  attack  required  by  the  pro- 
gram in  cases  considered  thus  far  have  not  exceeded  the 
-28°  to  +28°  range,  and  no  attempt  has  been  made  to  extend 
the  tables  set  up  for  input  at  this  time. 

A flap  is  used  to  control  the  twist  on  the  CTR,  so  the  aero- 
dynamic coefficients  for  this  rotor  vary  with  three  factors: 
angle  of  attack,  Mach  number,  and  flap  deflection.  Tables 
were  compiled  and  presented  in  reference  17  for  c^,  cd  and  cm 

for  the  modified  NACA  23012  airfoil  with  and  without  a flap. 
The  tables  covered  angles  of  attack  from  0°  to  360°,  Mach 
numbers  from  0.3  to  0.8,  and  flap  deflections  from  -10°  to 
+10°.  For  the  sections  of  the  blade  without  a flap,  the  co- 
efficients will  depend  on  angle  of  attack  and  Mach  number. 
Equation  (16)  can  easily  be  extended  to  include  an  interpo- 
lation on  Mach  number,  so 


I 

F 


' • 1 3a  1 '4g  ■ 


ct(a,M)  = (l-da)  cl(ai,Mj)  + Vt(ai+l'Mj)  + dM 

[<l-<Ja>  ct(a.,Mj+1)  + 4aOl<<-i+l'Mj+l)  -(1'da> 


or 


ct(a,M)  = (l-da)  (l-drf)  ci(ai.Mj)  +^(1-^)  ^(a^M.) 


+ (l-do)  dMC£(ai/Mj+i'  + d0dMCfc {ai+l'Mj+l} 


(17) 


where 


M-M. 


d..  = 


M Mj+l‘Mj 


For  Mach  numbers  below  0.3,  the  flow  is  incompressible,  so 
the  value  given  for  M = 0.3  was  used  at  all  Mach  numbers  be- 
low 0.3.  Values  of  the  coefficients  for  Mach  numbers  above 
0.8  were  obtained  by  extrapolation  from  the  value  given  at 
M = 0.8. 


Finally,  extending  Equation  (17)  to  include  interpolations 
for  flap  deflection  and  using  a subscript  "i,j,k"  to  indicate 
a tabulated  value  at  a.,  M.,  and  6,  will  yield 

1 1 rv 


Cfc (° ,M, fij)  * 


(l-da) (l-dM) (i-djlc,  + da(l-dM) (l-df)c 

»k  i + 1, j ,k 


+ (1-d«>dM(1-df,C.i  1+1  k + »V«.  ■ 

i,3+l,k  :,3,k+l 


+ ',adM(1‘df)ct  + da(1-Vdfci 

£ 4i+l,j+l,k  0 K f li+l,j,k+l 


+ (1-d  )d  dfc  + d dMd,c, 

a M £ ‘i,j+i,k+i  ° M 1 *i+i,j+i.k*i 


(18) 
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^ ‘i'-iERW1^  5UW*?.Wl! 


^T^rrr 


where 


6f-6f, 


^f  ~^f 
rk+l  rk 


Lumped  Loads 


Conversion  of  the  aerodynamic  loads  to  a form  suitable  for 
response  calculations  is  done  by  lumping  the  distributed 
loads  at  the  mass  points  of  the  lumped  parameter  blade  model. 
The  distributed  aerodynamic  moment  is  transferred  to  quarter- 
chord;  then  the  distributed  loads  are  integrated,  using  a 
straight-line  approximation  between  load  points,  to  obtain, 
lumped  loads  at  the  mass  points.  The  drag  force  at  the  i-— 
mass  point,  for  example,  is 


i+1 


F'(r)  dr 


R. 


where  R^  and  are  midway  between  masses  pk  and  iik_^  and 

between  m^+^  and  nu,  respectively,  and  F^(r)  is  the  straight- 

line  approximation  to  the  distributed  drag  load.  Similarly 
the  lumped  moment,  mo  , and  lift  force,  f„  , are  computed 


zi 


from  Mq  and  F . Coordinate  transformations  are  then  applied 


which  result  in  loads  in  the  local  blade  coordinate  system, 
and  are  given  by 


vi 


* “ f COS  t.  - f sin  <&. 


wi 


F sin  <t> . - f cos  4>. 
zi  1 xi  1 


M, 


mo  + z f 
i ai  vi 
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t.,.. 


— ~1'"' ViMiMBiiiiil jjiajjgj 


a 


K%s& 


where  z = distance  of  the  elastic  axis  forward  of  quar- 
ai  ter  chord,  feet 

4.  = angle  between  chord  and  the  plane  normal  to  the 
shaft,  positive  for  airfoil  rotated  nose  up, 
radians 


The  force  at  the  i—  mass  station  and  azimuth,  in  the 
direction  of  blade  flapping  motions,  denoted  by  Q , is 

vij 

given  by 

2 

Q = m. fi  sin4.  . (h.  + e.cos  4.  .)  + m . 0 [e.  + 

1 XJ  1 X XJ  X Cj  X 

hicos(4>i_.  - e0j)/cos  *i  - 2mifi^i;.  eisin  4^  + Pv 


+ m.f  (cos  4.. cos  a - sin  4.. sin  a siniji..) 
i g ns  ns  rn 


Here,  the  first  term  is  the  steady  acceleration  force  on  the 

i — mass  due  to  rotation,  the  second  term  accounts  for  con- 
trol angle  acceleration,  the  third  term  is  a gyroscopic 

coupling  term,  the  fourth  term  is  the  aerodynamic  load  on  the 
fch  i.1. 

i — mass,  and  the  fifth  term  is  the  weight  of  the  i—  blade 
segment.  The  in-plane  moment  is  given  by 


Qv  * m.fl^  e.r.  - 2fl0  [I j sin  4.  . + m.e.h.  sin  80/cos  4..] 

X XX  v«  XJ  XXX  Xj 

• • • 

- 2f24> . . I o + 2m.  SU  . (v. .sin  4..  + w.  .cos  4,.) 

XJ  X X XJ  XJ  XJ  XJ 


and  the  in-plane  force  is 


2w. . = min2  cos  $ij(hi  + ei  cos  *ij>  “ mihiec  sin($ij  ' 
/cos  4^ j 

- 2m.SU. ¥ . .cos  4.  . + P - m.f  (sin  4.  .cos  a + 
i i 13  i]  w^j  1 g 13  s 

cos4 . .sin  a sin  1(1 .. ) 

X J 5 X J 
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Blade  Responses 

Generalized  Forces  and  Coordinates 


The  generalized  forces  acting  on  each  normal  mode  are 
computed  for  each  azimuth , according  to 


F (t)  = la  ? + I 

K JKJJ  i 


0 A (K)  f Q A (K) 
vi  vi  wi  wi 


^ Q*  A,  (K)  + Qf  A^  <K>  + Q Aq  <K> 
*i  *i  *i  h *i  [ 


where  c is  from  the  previous  iteration, 
J 


°kj  - " cd;  \ W \ <J)  £or  3 * x 


°KK  = 2 °K  UK 


av  s the  average  aerodynamic  damping  coefficient  for 

«■  0.1* 

the  K—  mode  (read  as  input) 


(K) 

A = the  mode  shape  quantity  for  the  "q"  type  of 
qi 


th 


clastic  deformation,  at  the  i — radial  location,  for 


the  K—  mode 


cn.  = the  torsional  damping  associated  with  motion 

ue 


(K) 


defined  by  Aa  ' . The  aRJ  terms  may  be  thought  of  as 


4> 


damping  coupling  terms. 


.th 


The  governing  equation  for  the  K-1-  generalized  coor- 
dinate, ?K,  is  given  by 


?K  + 2aKaiK  ?K  + UK  ?K  " PK(t) 


(19) 
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The  damping  term,  oR,  is  defined  in  terms  of  the  mode 


shape  quantities  and  as 


°K  * aK  + °D 


r*.  <k'12A 
§ *,  / 


It  should  be  noted  that  viscous-type  damping  of  any 
motion  (torsional,  flatwise,  or  ohordwise)  may  be  repre- 
sented by  a similar  expression.  The  solution  of  this 
equation  is  obtained  in  integral  form.  The  solution 
assumes  periodicity  of  both  forcing  function  and  response. 
The  integral  relations  for  the  generalized  coordinate, 

CR,  and  its  time  derivatives  are  given  in  references 

7 and  22. 

The  average  aerodynamic  damping  coefficient  term,  2aKu>K, 

occurs  on  both  sides  of  the  governing  equation  for  cR(t). 

Since  an  iterative  solution  method  is  used  in  computing 
a compatible  set  of  loads  and  responses,  this  term  ef- 
fectively cancels  at  convergence. 

Blade  Response  Quantities 
Response  variables  are  computed  from 

v(t)  -l  Av<K)iK(t) 


mi  - j VK)=K<t> 


and  similarly  for  w,  and 

where  v = flatwise  blade  deflection 
w = chordwise  blade  deflection 
4>  * torsional  blade  deflection,  radians 
41  = 3v/8r,  radians 
V = 3w/8r,  radians 


for  the  torsional  moment,  T,  and  similarly  for  other 
moments,  shears,  and  motions. 

j * DETERMINATION  OF  PERFORMANCE  PARAMETERS  AND  STABILITY 
| AND  CONTROL  DERIVATIVES 

! 

Transfer  of  Forces  and  Moments  to  the  Hub 

i A single  blade  analysis  is  conducted  by  the  program  during 

its  loads-response  iterations.  Upon  converging,  the  normal 
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and  chordwise  shears,  moments,  displacements,  velocities, 
etc.,  are  computed  for  each  mass  point  along  the  blade  and 
for  each  azimuthal  position.  A harmonic  analysis  of  the 
shears  and  moments  at  the  root  of  the  blade  is  then  carried 
out.  For  articulated  blades  the  blade  root  is  taken  to  be 
the  location  of  the  flapping  hinge (s),  and  for  cantilevered 
blades  (such  as  those  for  the  ABC  rotor)  the  root  is  the  hub. 
The  axis  system  for  the  shears  and  moments  is  normal  to  and 
in  line  with  the  chord  of  the  airfoil,  as  shown  in  Figure  7. 
Resolving  these  shears  and  moments  into  a system  in  line 
with  and  normal  to  the  shaft  then  requires  that 


M 

z 

= Mz 

cos 

e 

+ 

M 

y 

sin 

0 

M 

= M 

cos 

e 

sin 

0 

y 

y 

z 

>i 

> 

y 

cos 

e 

+ 

V 

z 

sin 

0 

V 

= V 

cos 

e 

V 

sin 

0 

z 

z 

y 

A Fourier  series  representation  for  the  hub  shears  and  mo- 
ments for  each  individual  blade  can  be  written  as 


M = M + 7 Iml  COS  n\l>  + M sin 

o n=l  - c_  s -1 

n ” 


n 


with  similar  relations  for  M , V , and  V . 

y y c 


Following  the  tables  and  methods  given  in  reference  25,  the 
Fourier  coefficients  for  the  root  shears  and  moments  of  each 
blade  are  used  to  determine  the  combined  effect  of  all  of  the 
blades  on  the  forces  and  moments  transmitted  to  the  hub. 

Thus,  we  find  that 


thrust  = N,  V 

b y. 


drag  = -1/2  Nfe  Vz 


pitching  moment  = 1/2  N 


rolling  moment  = 1/2 


N Jk 


■»[\ + V \ XH 


- M, 


- V 


' 


sign(fi) 


The  rotor  torque  is  determined  from  distribution  of  aerody- 
namic drag  along  the  radius  of  the  blade,  so 


R 

torque  = N^J  d(r)  r dr 


V1 

- l/6Nb  l 

j-0 


d5(rj+i+2rj>  + dj+i(2rj+itrj>  (rj+rrj> 


where  is  the  number  of  blades,  d^  is  the  drag  acting  at  r^ 
and  N is  the  number  of  radial  load  points. 


Stability  and  Control  Derivatives 

For  stability  and  control,  the  C.G.  of  the  aircraft  is  the 
reference  point  to  be  used.  Thus  the  forces  and  moments  at 
the  rotor  hub(s)  must  be  transferred  to  the  C.G.  As  one 
example,  this  transfer  is  especially  important  for  pitch 
control  for  a tandem-rotor  system  where  the  pitching  moment 
consists  almost  entirely  of  the  thrust  of  each  rotor  times 
its  offset  from  the  C.G.  Thus,  pitch  control  for  a tandem 
rotor  system  is  achieved  through  the  differential  collec- 
tive pitch  between  the  fore  and  aft  rotors. 

The  stability  and  control  derivatives  are  computed  by  per- 
turbing the  shaft  angle  or  a blade  pitch  setting  from  its 
trim  value  and  then  dividing  the  resulting  change  in  the 
performance  parameters  by  this  perturbation.  For  example, 
consider  the  principal  thrust  control  for  a single-rotor 
helicopter.  In  this  case  the  collective  pitch  setting  would 
be  changed  from 


0 to  0 + 60  , 

o o o 

and  the  thrust  (as  well  as  all  other  forces  and  moments)  is 
recomputed.  Then  the  thrust  control  is  given  by 

T (0  + 60  ) -T ( 0 ) 

3T  _ o o o 

3 0 60 

0 0 

The  perturbation  variables  (e.g.,  60  , 6Alf  6B^,)  must  be 

large  enough  so  that  round-offs  won't  cause  an  error  by  try- 
ing to  find  a small  difference  between  two  large  numbers  and 


v ' iwjiyju. 


y«t  not  so  large  that  nonlinearities  will  be  introduced. 
Perturbations  of  0.2  deg  to  0.5  deg  have  been  found  to  fit 
the  forementioned  requirements  well.  Changing  the  lateral 
and  longitudinal  cyclic  pitch  settings  provides  the  neces- 
sary values  for  finding  the  roll  and  pitch  control  which  is 
available  from  the  rotor  on  a single-rotor  system. 


For  dual-rotor  systems,  differential  pitch  settings  can  be 
used  between  the  rotors  in  order  to  achieve  the  desired  con- 
trol. For  example,  differential  collective  pitch  is  used 
for  pitch  control  on  a tandem  rotor.  Some  differential  col- 
lective may  be  required  to  achieve  trim  conditions. 

Changing  this  differential  pitch  then  provides  the  pitch  con- 
trol from 


3PM 
3 A6 


PM 


A0 


+ sue  ) 

0 


-PMUeQ) 


SUIT) 


Thus  it  is  apparent  that  in  order  to  find  the  stability  and 
control  derivatives  for  rotor  systems  in  general,  one  must 
(1)  determine  the  values  for  the  performance  parameters 
when  the  aircraft  is  trimmed;  (2)  perturb  the  shaft  angle, 
blade  pitch  settings,  and  differential  pitch  settings  (for 
dual-rotor  systems)/ one  at  a time, and  recompute  the  perform- 
ance parameters  corresponding  to  these  new  settings;  and 
then  (3)  divide  the  performance  changes  by  the  perturbations 
made  on  the  settings  to  find  the  various  derivatives.  The 
derivatives  of  all  of  the  performance  parameters  with  re- 
spect to  all  attitude  and  pitch  variables  will  then  have 
been  determined.  The  principal  control  derivatives  of  in- 
terest for  the  various  rotor  systems  are  given  in  the 
following  table. 
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Rotor  type 

thrust 

pitching 

moment 

rolling 

moment 

yaw 

moment 

tingle 

e 

0 

Bi 

Ai 

* 

T 

TR 

tandem 

0 

0 

A0 

0 

Ai 

AA! 

toaxial 

0 

Ai 

ABj 

A0 

tounter  rotating 

| 0 

0 

tide-by-side 

0 

0 

Bi 

A0 

0 

ABj 

is  tail  rotor  thrust 


RESULTS  AND  DISCUSSION 


The  cases  considered  in  this  investigation  included  (1)  re- 
presentative flight  conditions  for  five  different  types  of 
rotor  systems  operating  at  advance  ratios  from  0.0  to  0.5, 

(2)  three  cases  for  the  ABC  rotor  at  p = 0.208  for  compari- 
son with  experimental  results,  and  (3)  three  cases  in  hover 
for  the  ABC  rotor  system.  Representative  characteristics 
of  the  five  rotor  systems  are  given  in  Tables  1-4,  and  the 
flight  conditions  ..nd  trim  settings  for  each  case  are  tabu- 
lated in  Table  5.  Discussion  of  the  results  obtained  for 
the  various  rotor  systems  is  separated  into  three  areas: 

1)  The  effect  of  the  deformed  wake  on  the  performance 
parameters  for  the  five  rotor  systems  at  the  par- 
ticular flight  condition  considered  for  each. 

2)  The  effect  of  the  deformed  wake  on  the  performance 
of  the  ABC  rotor  system  at  p = 0.208  and  in  hover. 

3)  The  effect  of  the  deformed  wake  on  the  stability 
and  control  derivatives 

a)  for  the  conventional,  tandem,  TRAC,  and  CTR 
rotor  systems  at  their  particular  flight 
conditions 


b)  for  the  ABC  rotor  system  in  forward  flight 
(including  comparisons  with  experimental  re- 
sults) and  in  hover. 

EFFECTS  OF  DEFORMED  WAKE  ON  THE  PERFORMANCE  CHARACTERISTICS 


OF  REPRESENTATIVE  ROTOR  SYSTEM 


The  five  types  of  rotor  systems  analyzed  in  this  study  were 
the  conventional,  single  rotor,  tandem,  TRAC,  CTR,  and  ABC 
systems.  The  characteristics  of  each  representative  rotor 
system  are  given  in  Tables  1 through  4.  The  conver  ional 
rotor  was  the  articulated  rotor  system  of  the  H-34  helicop- 
ter as  described  in  reference  26.  The  tandem-rotor  system 
used  two  H-34  rotors  orientated  to  provide  an  80%  overlap 
between  the  forward  and  aft  rotors  with  a 5.67-ft  (0.20R) 
vertical  separation.  The  representative  Telescoping  Rotor 
Aircraft  system  was  simulated  by  using  the  XH-51  helicopter 
ElacTes  with  a 40%  cutout.  The  characteristics  of  the  Con- 
trollable Twist  Rotor  were  obtained  from  reference  17.""  The 
aerodynamic  center  was  assumed  to  lie  at  25%  of  the  total 
chord  and  thus  was  2 inches  behind  the  blade  elastic  axis 
for  the  blade  sections  with  the  flap.  In  addition,  it  was 
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assumed  that  the  first  torsional  frequency  of  the  blade  was 
4/rev.  A cambered  airfoil  is  used  on  the  CTR,  so  the  zero- 
lift  angle  was  non-zero.  This  zero-lift  angle  and  its 
variation  with  flap  deflectio.i  are  given  in  Table  4 along 
with  the  spanwise  location  and  chord  length  of  the  flap. 

The  Advancing  Blade  Concept  system  consists  of  two  counter- 
rotating,  coaxial,  cantilevered  rotors.  The  blades  which 
were  built  for  wind  tunnel  testing  are  described  in  refer- 
ence 16.  The  nonlinear  built-in  twist  and  the  radial 
distribution  of  the  airfoil  shape  for  the  ABC  are  given  in 
Tables  2 and  3. 

The  flight  conditions  shown  in  Table  5 were  chosen  so  as  to 
provide  a representative  case  for  each  rotor  system  and  to 
cover  advance  ratios  from  hover  to  p = 0.5.  The  trim 
settings  for  blade  pitch  and  the  blade  flapping  response 
given  in  Table  5 were  obtained  from  reference  26  for  the 
H-34  helicopter,  from  trim  runs  with  the  U.  S.  Army's 
Rotorcraft  Flight  Simulation  Program,  C-81,  for  the  TRAC 
and  tandem  systems,  and  from  reference  17  for  the  CTR. 
Reference  16  provided  blade  pitch  settings  which  were  the 
average  value  between  the  upper  and  lower  rotors  for  the  ABC. 
These  average  values  were  used  as  starting  points  for  a manual 
trial  and  error  trim  procedure  that  was  followed  to  find  the 
individual  pitch  settings  of  each  rotor. 

The  performance  parameters  calculated  by  the  program  de- 
veloped herein  with  and  without  deformed-wake  effects 
included  are  given  in  Table  6 for  the  conventional,  TRAC, 
tandem,  and  CTR  rotor  systems. 

A uniform  induced  velocity  distribution  was  used  for  cases 
run  without  deformed-wake  effects  included.  This  uniform 
inflow  is  computed  from  momentum  considerations,  and  no 
wake  effects  are  included  in  this  type  of  model.  When  de- 
formed-wake effects  were  included,  the  radial  and  azimuthal 
distribution  of  the  induced  velocity  was  computed  from  the 
known  positions  and  strengths  of  the  vortices  in  the  wake. 
Wake  effects,  then,  show  up  when  the  results  obtained  with 
the  nonuniform  wake-induced  inflow  are  compared  with  the 
results  gained  assuming  uniform  (i.e.,  momentum)  inflow. 

The  weight  of  the  H-34  during  the  flight  tests  varied  be- 
tween 11,200  and  11,805  lb  (see  Ref.  26).  The  thrust  of 
11,878  lb  shown  in  Table  6 with  wake  effects  included  and 
using  measured  trim  control  settings  is  in  excellent 
agreement  with  the  measured  values.  The  wake  did  not  have 
a significant  effect  on  any  of  the  performance  character- 
istics for  this  particular  flight  condition  for  the  H-34 


since  it  represents  operations  at  moderate  values  for  the 
loading,  advance  ratio,  and  shaft  tilt  of  a single-rotor 
system  where  the  wake  effects  would  be  expected  to  be  small. 

Trim  settings  for  the  TRAC  and  tandem-rotor  systems  were 
obtained  by  running  the  trim  option  of  the  C-81  program  for 
each  system.  The  performance  parameters  computed  by  the 
program  developed  herein  using  uniform  inflow  are  in  fair 
agreement  with  those  computed  by  the  C-81  program.  Better 
agreement  may  have  resulted  if  blade  mode  shapes  had  been 
input  to  the  C-81  program  to  account  for  blade  flexibility 
during  its  trim  iterations. 

The  wake  remained  close  to  the  rotor  in  the  hover  condition 
used  for  the  TRAC,  and  results  shown  in  Figure  8 indicate 
a drop  of  8.7%  in  thrust  and  24.8%  in  the  power  required 
between  the  case  for  uniform  inflow  and  the  case  using  the 
predicted  distribution  of  the  wake-induced  velocity  for  the 
same  control  settings.  Retrimming  the  rotor  with  wake  ef- 
fects included  would  require  that  the  collective  pitch  be 
increased  in  order  to  recover  the  thrust  lost  due  to  wake 
effects.  The  resulting  increase  in  required  power  might 
bring  it  back  up  to  the  1029  hp  which  was  needed  for  uniform 
inflow. 

The  trust  produced  and  power  required  for  the  CTR  with  and 
without  wake  effects  included  are  presented  in  Figure  9. 
Although  the  thrust  decreased  by  7.1%  with  wake  effects  in- 
cluded, the  power  requirements  increased  by  15.4%.  More 
flight  conditions  with  subsequent  retrimming  with  wake  ef- 
fects included  need  to  be  considered  in  order  to  further 
investigate  and  explain  this  power  increase  with  thrust  loss 
when  wake  effects  are  included.  These  changes  are  quite 
significant  even  though  the  CTR  is  a single,  articulated 
rotor  system  and  the  advance  ratio  and  shaft  angle  are 
fairly  large,  seemingly  reducing  the  wake  induced  effects. 
This  apparent  reversal  is  probably  due  to  large  radial  and 
azimuthal  variations  in  loading  for  the  CTR,  especially 
along  the  flap.  This  variation  would  produce  shed  and 
trailing  vortices  for  the  CTR  which  are  quite  strong  and 
thus  cause  strong  wake  effects. 

The  overlap  of  the  rotors  on  the  tandem  system  produced  a 
strong  wake  interaction  over  the  aft  portion  of  the  lower 
rotor  and  over  the  forward  portion  of  the  aft  rotor.  The 
wake  from  the  aft  rotor  swept  through  the  aft  portion  of 
the  forward  rotor  and  significantly  reduced  the  induced 
downwash  in  this  area.  The  result  was  a 21.9%  increase  in 
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thrust  and  42.5%  increase  in  power  required  for  the  forward 
rotor  when  compared  to  results  obtained  assuming  uniform  in- 
flow. The  wake  from  the  forward  rotor  moved  under  the  aft 
rotor  and  increased  the  downwash  for  the  aft  rotor.  The 
wake  interaction  effect  was  smaller  on  the  aft  rotor  than  it 
was  on  the  forward  rotor  since  the  wake  is  farther  away  from 
the  blades,  and  therefore  th2  thrust  and  power  drops  were 
only  8.6%  and  11.9%,  respectively.  These  thrust  and  power 
changes  for  both  rotors  as  well  as  the  pitching  and  rolling 
moments  at  the  hub  of  each  rotor  are  illustrated  in  Figure 
10. 


EFFECTS  OF  THE  DEFORMED  WAKE  ON  THE  PERFORMANCE  CHARACTERIS- 
TICS OF  THE  ABC  ROTOR  CONFIGURATION 

The  performance  parameters  which  were  computed  for  each 
rotor  in  the  ABC  system  are  listed  in  Table  7 for  the  vari- 
ous flight  conditions.  The  results  for  the  forward  flight 
cases  are  shown  in  Figures  11  through  14,  and  the  hover 
cases  are  shown  in  Figures  16,  17,  and  18.  Since  the  ABC 
is  a coaxial  rotor  system,  the  loading  distribution  was 
significantly  different  on  the  lower  rotor  when  comparing 
the  results  obtained  with  and  without  wake  effects  included. 
In  hover,  the  wake  from  the  upper  rotor  passes  through  all 
of  the  lower  rotor  and  creates  a highly  nonuniform  upwash- 
downwash  distribution.  The  tunnel  tests  in  forward  flight 
were  conducted  mainly  with  the  shaft  being  vertical  or 
tilted  aft.  This  caused  the  wake  from  the  upper  rotor  to 
remain  primarily  above  the  lower  rotor,  and  in  some  condi- 
tions the  wake  from  the  lower  rotor  also  moved  up.  In  these 
flight  conditions,  a large  thrust  increase  would  be  obtained 
on  the  lower  rotor,  and  it  would  go  into  an  autorotation 
mode  and  produce  power  rather  than  use  it. 

In  Case  1,  the  rotors  were  operating  at  an  advance  ratio  of 
0.466  and  the  shaft  was  vertical.  No  significant  change 
occurred  in  the  thrust  or  power  of  the  upper  rotor  when  the 
effects  of  the  deformed  wake  were  included.  The  changes 
shown  in  Figure  II  for  the  lower  rotor  were  a 24.5%  in- 
crease in  thrust  and  a 40.4%  decrease  in  power.  The 
pitching  and  rolling  moments  for  both  rotors  increased,  and 
the  net  result  was  a change  of  the  total  pitching  moment 
from  zero  to  +5100  ft-lb  and  a change  of  the  total  rolling 
moment  from  -4000  to  +4800  ft-lb. 

The  advance  ratio  for  Cases  2,  3,  and  4 for  the  ABC  was 
0.208,  and  the  shaft  was  tilted  aft  by  4°  for  Cases  2 and 
3 and  8°  for  Case  4.  The  result  of  the  calculations  con- 
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ducted  for  these  flight  conditions  indicated  a significant 
increase  in  thrust,  especially  for  the  lower  rotor  when  the 
wake  effects  were  included  For  the  lower  rotor,  the 
change  from  power  required  (i.e.,  positive  horsepower)  with 
uniform  downwash  to  power  produced  (i.e.,  negative  horse- 
power) with  the  actual  wake-induced  upwash  distribution  is 
obvious  in  Figures  12,  13,  and  14.  This  power  switch  is 
due  to  the  aft  shaft  tilt  which  causes  the  wakes  from  both 
rotors  to  be  swept  above  the  lower  rotor.  This  causes  an 
upwash  on  the  lower  rotor,  thus  placing  it  in  an  autorota- 
tive  mode  of  operation. 
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In  Case  2 the  shaft  angle  was  4°  aft  and  the  rotors  were 
trimmed  to  yield  a total  thrust  of  14,700  lb.  This  smaller 
shaft  angle  and  thrust  caused  the  wake  to  remain  nearby, 
so  it  had  a significant  effect  on  both  rotors.  The  results 
shown  in  Figure  12  for  Case  2 indicate  an  increase  in 
thrust  of  15.0%  for  the  upper  rotor  and  44.2%  for  the  lower 
rotor.  The  power  required  by  the  upper  rotor  decreased  by 
54.1%  when  the  wake  effects  were  included,  and  the  power 
for  the  lower  rotor  changed  sign  (i.e.,  negative  horsepower). 
The  pitching  and  rolling  moments  increased  for  both  rotors 
with  wake  included.  The  net  result  was  a change  from 
-2300  ft-lb  to  +3350  ft-lb  for  the  pitching  moment.  The 
net  rolling  moment  changed  from  -7400  ft-lb  to  +7700  ft-lb. 


The  thrust  was  increased  to  22,000  lb  for  Case  3.  This 
caused  the  wake  to  move  farther  away  from  the  upper  rotor, 
and  thus,  except  for  the  rolling  moment,  the  wake  had  only 
a slight  effect  on  the  performance  parameters  for  the 
upper  rotor.  The  wake  effects  for  Case  3 caused  an  up- 
wash on  the  lower  rotor,  and  the  result,  as  illustrated  in 
Figure  13,  was  a 29.2%  increase  in  thrust  for  the  lower 
rotor  and  a change  in  sign  for  the  power  required.  The 
thrust  and  power  remained  essentially  unchanged  for  the 
upper  rotor.  The  left  roll  produced  by  the  upper  rotor  de- 
creased by  22.9%,  and  the  right  roll  produced  by  the  lower 
rotor  increased  by  33.8%.  The  net  change  in  rolling  moment 
then  was  from  -11,250  ft-lb  to  +12,000  ft-lb. 

The  shaft  angle  was  increased  from  4°  to  8°  in  moving  from 
Case  3 to  Case  4.  Again  for  the  upper  rotor,  only  the 
rolling  moment  changed  significantly.  An  upwash  occurred 
on  the  lower  rotor;  as  a consequence,  the  results  given  in 
Figure  14  for  Case  4 show  a 35.4%  increase  in  thrust  and  a 
change  from  130  horsepower  to  -500  horsepower  for  the 
lower  rotor.  A very  large  increase  also  occurred  in  the 
pitching  moment  produced  by  the  lower  rotor.  The  left  roll 
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of  the  upper  rotor  decreased  by  13.5%,  and  the  right  roll 
of  the  lower  rotor  increased  by  37.4%.  The  net  result  was 
a change  from  8100  ft-lb  oi  left  roll  with  uniform  inflow 
to  8100  ft-lb  of  right  roll  with  wake  effects  included. 

The  results  for  the  performance  parameters  show  that  the 
rotor  system  should  be  retrimmed  when  wake  effects  are  in- 
cluded. The  net  effect  due  to  inclusion  of  the  wake  was 
generally  a large  increase  in  the  net  thrust,  pitching 
moment,  and  rolling  moment.  The  blade  pitch  settings  for 
both  rotors  need  to  be  changed  to  bring  these  parameters 
back  to  the  desired  trim  values.  With  the  lower  rotor  in 
autorotation,  the  torques  from  the  rotors  will  join  instead 
of  cancel  each  other.  As  a consequence,  a significant 
amount  of  differential  collective  pitch  may  be  required  in 
order  to  restore  yaw  equilibrium. 

The  wake  geometry  and  the  wake  influence  coefficients  for 
Cases  2 and  4 were  used  with  the  flight  conditions  and  trim 
settings  of  Case  3 in  an  effort  to  ascertain  the  sensitivity 
of  the  calculations  to  variations  in  the  wake  geometry. 

Case  2 had  required  a total  thrust  of  14,680  lb  instead  of 
the  21,980  lb  for  Case  3.  The  shaft  angle  was  changed  from- 
4°  to  8°  between  cases  3 and  4,  respectively.  The  differ- 
ences in  the  wake  geometries  caused  by  these  differences 
in  the  flight  conditions  did  not  produce  any  significant 
changes  in  the  performance  parameters,  except  perhaps  in 
the  power  requirements  for  the  lower  rotor.  Tne  horse- 
power ranged  from  -225  horsepower  to  -280  horsepower  among 
the  cases  for  the  different  wakes,  and  this  represents  a 
variation  of  +12.2%  from  the  average  value  for  the  three 
cases.  The  thrust  and  rolling  moment  variations  shown  in 
Figure  15  are  all  less  than  2%  for  the  different  geometries. 

Thrusts  of  10,000,  15,000,  and  20,000  pounds  with  zero  shaft 
tilt  angle  were  considered  for  the  ABC  in  hover.  The  re- 
sults for  the  thrust  and  power  of  each  rotor  are  shown  in 
Figures  16,  17,  and  18.  Since  the  induced  effect  of  the 
lower  rotor  moved  the  upper  rotor  vortices  inward,  the  wake- 
induced  velocity  showed  upwash  through  the  outer  portions 
of  the  lower  rotor  disc  Ln  these  hover  cases.  As  a result, 
the  thrust  for  the  i’v'-v  -otor  increased  by  31%  to  38%  for 
these  various  cases.  Th.-  ihrast  increase  for  the  upper 
rotor  was  between  I.'*  and  . since  the  wake-induced  down- 
wash  was  generally  small-,  r chan  the  value  used  for  uniform 
inflow.  The  effect  of  the  wake  on  the  power  required  by 
the  upper  rotor  was  very  small  in  all  cases,  but  the  effect 
on  the  power  required  by  the  lower  rotor  was  quite  large. 
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Again,  because  of  upwash  through  part  of  the  lower  rotor,  the 
power  required  by  the  lower  rotor  decreased  by  about  500 
horsepower  in  each  case.  This  meant  a sign  change  for  the 
horsepower  for  T = 10,000  lb,  a 66.9%  decrease  in  horsepower 
for  T = 15,000  lb  and  a 30.4%  decrease  for  T = 20,000  lb. 

The  variation  of  the  thrust  from  either  rotor  and  the  power 
of  the  upper  rotor  was  less  than  3%  when  uhe  wake  geometries 
for  T = 10,000  and  T = 20,000  lb  were  used  with  the  trim 
collective  and  cyclic  pitch  settings  for  T = 15,000  lb.  The 
range  of  power  required  by  the  lower  rotor  was  +95  horsepower 
(i.e.,  51%)  from  the  average  for  the  three  cases  shown  in 
Figure  19  with  wake  effects  included. 

Comparing  the  results  obtained  for  the  performance  parameters 
with  and  without  wake  effects  included  for  the  ABC  has  shown 
that  use  of  the  correct  wake-induced  velocity  distribution 
is  very  important  for  coaxial  counterrotating  rotor  systems. 
This  is  especially  true  for  determining  the  forces  and  mo- 
ments generated  by  the  lower  rotor.  The  wake  effects 
generally  caused  an  autorotative  mode  of  operation  for  the 
lower  rotor,  thus  increasing  its  thrust  and  decreasing,  or 
even  reversing,  the  power  required  by  the  lower  rotor.  For 
the  forward  flight  cases  considered  here,  a large  right 
rolling  moment  was  introduced  when  wake  effects  were  in- 
cluded, The  pitching  moment  was  generally  more  nose-up  with 
wake  effects  included.  Changing  the  blade  root  pitch  set- 
tings in  a retrim  procedure  would  be  required  in  order  to 
bring  the  thrust  and  moments  back  to  the  desired  trim 
values  with  wake  effects  included. 


WAKE  EFFECTS  ON  STABILITY  AND  CONTROL  DERIVATIVES 

Conventional,  TRAC,  CTR,  and  Tandem  Rotor  Systems 

The  wake  had  only  a slight  effect  on  most  of  the  stability 
and  control  derivatives  for  the  articulated,  single-rotor 
systems  for  the  flight  conditions  considered  in  this  study. 

For  a moderate  advance  ratio  and  shaft  tilt  angle,  the  wake  was 
quickly  carried  away  from  the  rotor  due  to  the  free  stream 
velocity  components  normal  to  and  parallel  with  the  tip  path 
plane.  As  a consequence,  the  effect  of  uniform  or  wake- 
induced  velocity  distributions  on  the  angle  of  attack 
becomes  overshadowed  by  the  effect  of  the  throughflow  velo- 
city, V sin  ag. 

The  values  computed  for  these  single-rotor  systems  with  and 
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without  wake  effects  included  are  given  in  Table  8.  The  de- 
rivative of  thrust  with  respect  to  collective  pitch  angle 
OT/36q)  increased  from  187  lb/deg  to  365  lb/deg  for  the 

H-34  and  from  149  lb/deg  to  290  lb/deg  for  the  TRAC  when 
deformed -wake  effects  were  included.  The  change  in  yawing 
moment  (i.e.,  rotor  torque)  due  to  collective  pitch  for  the 
TRAC  was  decreased  by  20.3%,  but  this  derivative  is  a coup- 
ling term  for  single-rotor  systems  since  the  principal  yaw 
control  would  come  from  the  tail  rotor  thrust. 

The  strong  wake  interaction  which  occurs  for  the  tandem- 
rotor  system  caused  a significant  change  in  the  stability 
and  control  derivatives  for  the  tandem  rotor.  The  values 
for  the  more  important  derivatives  for  the  tandem  rotor  with 
and  without  wake  effects  included  are  given  in  Table  9. 

The  thrust  and  pitching  moment  (about  the  C.G.)  controls 
shown  in  Figure  20  were  reduced  by  34.8%  and  25.9%,  respec- 
tively, when  the  wake  effects  were  included.  With  the 
deformed-wake  included,  an  increase  in  the  thrust  produces 
a stronger  wake,  and  thus  a greater  downwash.  Some  collec- 
tive pitch  is  required  to  overcome  this  larger  downwash  and 
is  lost  as  far  as  thrust  control  is  concerned.  As  a conse- 
quence, the  thrust  control  derivative  with  respect  to  col- 
lective pitch  is  reduced  by  including  wake  effects.  This 
reduction  is  evident  in  Figure  20.  The  same  interdependence 
between  thrust  and  wake-induced  downwash  reduces  the  pitch 
control  since  the  pitching  moment  is  mainly  composed  of  the 
differential  thrust  between  the  two  rotors  times  their  hori- 
zontal offsets  from  the  C.G. 

ABC  Rotor  System,  Hover  and  Forward  Flight 

Advance  ratios  of  0.208  and  0.466  were  considered  for  the 
ABC  so  that  comparisons  could  be  made  with  the  experimental 
values  given  in  reference  16  for  these  advance  ratios.  The 
trim  settings  and  the  forces  and  moments  generated  by  each 
individual  rotor  were  not  given  in  reference  16.  Since  the 
trim  settings  he  1 to  be  obtained  by  trial  and  error  using 
iterative  single  runs  of  the  blade  loads-response  program, 
only  three  cases  at  p = 0.  208  and  one  case  at  ;j  = 0.466 
could  be  completed  in  the  time  available  the  calcula- 
tions. These  cases  did  include  variatio  n shaft  angle, 

thrust,  and  flight  speed,  and  the  effect  these  param- 

eters on  the  stability  and  control  derivatives  will  be 
discussed  here. 

The  thrust  and  moments  produced  by  each  rotor  were  not  known 
The  total  thrust  is,  of  course,  the  sum  of  the  thrusts  of 
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the  upper  and  lower  rotors.  The  total  rolling  moment  is  com- 
posed of  the  thrust  from  each  rotor  times  its  lateral  lift 
offset.  The  lift  offset  given  in  reference  16  for  each 
flight  condition  was  assumed  to  apply  to  both  rotors.  Thus, 
the  thrust  from  each  rotor  was  determined  from  the  following 
two  equations: 

Tl  + Tu  T total 

rl(Tl  - T ) = RMtotal 


where  rL  is  the  lateral  lift  offset. 


Each  rotor  was  trimmed 


out  to  produce  zero  pitching  moment.  Considering  each  rotor 
separately  then,  the  thrust,  rolling  moment,  and  pitching 
moment  were  the  dependent  variables  which  were  controlled  by 
the  blade  pitch  settings  (collective,  lateral  cyclic,  and 
longitudinal  cyclic) . 

Since  an  automated  trim  routine  was  not  available  for  the 
ABC  rotor  system,  a manual  trim  procedure  was  employed  which 
involved  varying  the  blade  pitch  settings  and  plotting  the 
resulting  performance  parameters  versus  these  pitch  settings. 
Several  runs  are  required  in  order  to  obtain  the  points 
needed  for  these  plots.  It  is  fairly  difficult  to  obtain 
the  trim  settings  assuming  uniform  inflow,  and  the  procedure 
becomes  too  complex  when  wake  effects  are  included.  Lacking 
an  automatic  trim  procedure,  the  upper  and  lower  rotors  were 
therefore  not  in  trim  when  the  wake-induced  velocity  dis- 
tribution was  used  with  the  pitch  settings  calculated 
assuming  uniform  inflow.  Since  the  true  trim  forces  for  each 
rotor  were  not  known,  no  attempt  was  made  to  retrim  the 
rotors  by  hand  with  wake  effects  included. 

Hover  Cases 


Thrusts  of  10,000,  15,000,  and  20,000  pounds  were  considered 
in  hover.  The  values  for  the  control  derivatives  for  these 
cases  with  and  without  wake  effects  included  are  given  in 
Table  10.  The  effect  of  thrust  on  the  control  derivatives 
with  and  without  wake  effects  included  is  shown  in  Figure  21. 
The  values  for  each  derivative  have  been  normalized  by  using 
the  value  at  T = 15,000  lb  as  a base. 

The  same  trends  exist  with  respect  to  thrust  with  and  without 
wake  included  for  the  thrust  control  (3T/36  ).  The  thrust 
control,  however,  was  more  sensitive  to  thr&st  changes  with 
wake-induced  velocities  than  it  was  with  uniform  inflow. 
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Also,  reductions  of  up  to  20  percent  occurred  in  the  thrust 
control  when  the  uniform  inflow  was  changed  to  the  wake- 
induced  velocity  distribution. 

The  results  for  the  roll  control  ORM/aABj)  show  the  same 
trends  with  and  without  wake  effects  included.  In  either 
case,  greater  roll  control  is  possible  with  greater  thrust. 

In  changing  from  uniform  inflow  to  wake-induced  velocities, 
the  roll  control  was  reduced  by  18  to  20  percent.  The  pitch 
control  through  lateral  and  longitudinal  cyclic  pitch  de- 
creases slightly  as  the  thrust  is  increased  with  uniform  in- 
flow. With  wake  effects  included,  the  trends  with  respect 
to  thrust  are  mixed. 

Good  yaw  control  by  use  of  differential  collective  pitch 
(i.e.,  about  7000  ft-lb/deg)  was  computed  for  the  ABC  rotor 
system  in  hover.  However,  the  trend  with  increasing  thrust 
for  wake-induced  velocities  was  opposite  that  which  was  pre- 
dicted assuming  uniform  inflow.  With  uniform  inflow,  the 
yaw  control  increased  with  thrust,  but  it  decreased  with 
thrust  when  the  wake  effects  were  included.  As  was  noted 
earlier,  an  autorotation  condition  exists  for  the  lower 
rotor  when  the  actual  wake  effects  are  considered,  and  this 
effect  becomes  stronger  as  the  thrust  increases.  Thus  the 
yaw  control  trend  is  reversed  when  wake-induced  velocities 
are  used.  The  sensitivity  of  the  derivatives  to  perturba- 
tions on  the  wake  geometry,  as  shown  in  Figure  22,  was 
fairly  large  for  these  hover  cases.  Due  to  this  sensitivity, 
the  wake  geometry  should  be  recomputed  with  any  change  in 
the  flight  conditions  at  low  speeds  or  in  hover. 

Forward  Flight  Cases 

Under  an  Army  contract  with  Sikorsky  Aircraft,  the  full- 
scale  ABC  rotor  system  was  tested  in  the  NASA/Ames  Research 
Center  40-ft-x-80-ft  Wind  Tunnel.  Some  of  the  values  ob- 
tained for  the  stability  and  control  derivatives  from  these 
tests  as  well  as  those  computed  in  this  study  with  and  with- 
out wake  effects  included  are  given  in  Table  11.  The 
flight  conditions  which  changed  among  these  cases  were  shaft 
angle,  thrust,  and  flight  speed.  The  effect  of  these  flight 
variables  on  the  stability  and  control  derivatives  is  pre- 
sented in  Figures  23,  24,  and  25. 

1.  Effect  of  Shaft  Angle. 

For  cases  3 and  4,  the  advance  ratio  was  0.208, 

and  the  thrust  was  22,000  lb.  The  shaft  angle 
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changed  from  4°  to  8°  between  these  cases.  Divid- 
ing the  values  for  each  derivative  when  ag  = 8°  by 
its  value  when  og  = 4°  provided  the  normalized 

values  shown  in  Figure  23.  The  roll  and  pitch  con- 
trol derivatives  increased  by  10%  when  the  shaft 

angle  was  increased  from  4°  to  8°.  Increasing  the 
shaft  angle  created  a greater  autorotative  condi- 
tion, so  the  yaw  control  for  og  = 8°  was  half  (or 

less)  what  it  was  for  a =4°.  Using  differential 

s 

collective  pitch  as  yaw  control  in  forward  flight 
was  not  found  to  be  very  effective  either  in  the 
wind  tunnel  or  in  this  analysis.  The  cases  with 
uniform  inflow  did  not  follow  the  trends  shown  by 
the  experiments  for  shaft  angle  as  well  as  the 
cases  with  wake-induced  velocities  did. 

The  effect  of  the  shaft  angle  on  longitudinal  sta- 
bility is  destabilizing  in  all  of  the  forward  flight 
cases.  The  positive  values  obtained  for  3PM/3a 

s 

indicate  that  an  increase  in  shaft  angle  will  cause 
a nose-up  moment  and  thus  tend  to  drive  the  air- 
craft to  even  higher  shaft  angles.  Control  of  the 
pitching  moment  through  8Qor  Bj  was  virtually  un- 
changed assuming  uniform  inflow;  it  changed  by 
moderate  amounts  in  the  tunnel  tests,  and  the  cases 
with  wake-induced  velocities  showed  large  changes 
due  to  shaft  angle.  The  wake  location  relative  to 
the  rotors  is  sensitive  to  the  shaft  angle.  Signi- 
ficant changes  would  probably  occur  in  the  trim 
settings,  especially  for  the  lower  rotor,  if  the 
rotors  were  retrimmod  with  wake  effects  included. 
Thus,  better  agreement  might  be  obtained  between  the 
calculations  with  wake-induced  velocities  and  the  ex- 
perimental values  if  the  rotors  were  retrimmed. 

2.  Effect  of  Thrust 

The  thrust  for  Case  2 was  15,000  lb  and  the  thrust 
for  Case  3 was  22,000  lb.  The  gdvance  ratio  was 
0.208  and  the  shaft  angle  was  4°  in  both  Cases. 

The  values  for  the  derivatives  for  T = 22,000  lb 
were  normalized  by  their  value  at  T = 15,000  lb. 

These  results  are  shown  in  Figure  24.  The  roll  and 
pitch  control  derivatives  were  reduced  by  10  to  25 
percent  when  the  thrust  was  increased  from  15,000 
lb  to  22,000  lb.  The  prediction  of  the  thrust  effect 
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was  closer  to  the  experimental  results  with  wake- 
induced  velocities  used  than  with  uniform  inflow. 
Differential  collective  pitch  was  quite  ineffective 
in  yaw  control  for  T = 15,000  lb.  Increasing  the 
thrust  to  22,000  lb  more  than  doubled  the  torque 
required  by  each  rotor  and  then  differential  col- 
lective became  effective  as  a yaw  control  parameter. 
Thus,  the  large  change  shown  in  Figure  24  for 
3YM/3A6  is  due  in  part  to  the  low  base  values 
which  existed  at  T = 15,000  lb  for  the  experimental 
case,  the  uniform  inflow  case,  and  the  wake-induced 
velocities  case.  The  experimental  values  show  a 
35%  reduction  for  the  destabilizing  effect  of  shaft 
angle  and  about  a 20%  reduction  for  control  of 
pitching  moment  through  0q  or  Bj  due  to  the  in- 
creased thrust.  Little  change  occurred  for  any  of 
the  longitudinal  stability  derivatives  due  to 
thrust  when  uniform  inflow  was  used.  The  calcu- 
lations with  wake  effects  included  were  sensitive 
to  the  thrust,  but  the  thrust  computed  on  the  wake 
runs  was  quite  a bit  out  of  trim.  Retrimming 
these  cases  with  wake-induced  velocities  should 
bring  the  derivatives  into  better  agreement  with 
the  test  results. 


3.  Effect  of  Advance  Ratio. 


The  advance  ratio  was  changed  from  0.208  in  Case  2 
to  0.466  in  Case  1.  The  shaft  angle  also  changed 

from  4°  to  zero,  so  the  trends  to  be  noted  here 
are  not  purely  due  to  flight  speed.  The  thrust  in 
both  cases  was  15,000  lb.  The  values  for  the  de- 
rivatives at  m = 0.466  were  divided  by  their  value 
at  u = 0.208,  and  these  normalized  results  are 
shown  in  Figure  25. 

The  higher  speed  produced  much  greater  roll  con- 
trol since  the  difference  in  local  velocities  on 
the  advancing  and  retreating  sides  of  each  rotor 
will  add  to  the  difference  in  pitch  angle  (due  to 
Bj ) and  thus  increase  the  effectiveness  of  Bj. 

The  change  in  3RM/3ABj  due  to  advance  ratio  which 

was  computed  using  wake-induced  velocities  is 
closer  to  the  experimental  results  than  uniform 
inflow  was.  A 10%  increase  in  pitch  control  oc- 
curred with  the  increased  flight  velocity  for  both 
the  experiment  and  the  calculations  with  uniform 
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inflow.  Retrimming  the  cases  with  wake-induced 
velocities  should  bring  the  wake  results  into 
better  agreement  with  experiment. 

The  destabilizing  effect  of  shaft  angle  was  in- 
creased by  a factor  of  five  when  the  advance 
ratio  was  increased  from  the  value  in  Case  2 to 
that  in  Case  1.  Control  of  the  pitching  moment, 
either  through  0 or  B](  was  essentially  doubled 

with  this  increase  in  advance  ratio.  Good  agree- 
ment exists  between  the  computed  results  and  the 
experimental  values  for  all  of  the  derivatives 
except  c'PM/9 0 . increase  occurs  for  3PM/30q 

when  the  advance  ratio  increases.  In  Figure  25 
it  is  evident  that  a fair  to  poor  comparison 
exists  between  the  measured  change  in  3PM/30q  and 

the  values  computed  assuming  uniform  inflow.  A 
large  discrepancy  exists,  however,  when  the  wake- 
induced  velocity  distribution  was  used  since  the 
wake  geometry  calculations  are  sensitive  to  shaft 
angle  and  advance  ratio.  The  higher  velocity  and 
shallower  shaft  angle  in  Case  1 (as  compared  to 
Case  2)  made  the  wake  be  relatively  well  removed. 
As  a consequence,  the  change  in  3PM/30  due  to  a 

change  in  y from  0.208  to  0.466  was  overestimated 
by  the  computations  with  wake-induced  velocities. 
Retrimming  these  cases  using  wake-induced  velo- 
cities should  improve  the  results  obtained  with 
wake  effects  included. 

The  measurements  of  the  trends  which  correspond 
to  changes  in  shaft  angle,  thrust,  and  flight 
speed  using  the  programs  developed  for  this  study 
have  compared  fairly  well  with  the  experimental 
measurements.  More  detailed  comparisons  will  re- 
quire (1)  an  automated  trim  procedure  for  deter- 
mining trim  settings  using  wake-induced 
velocities  and  (2)  specific  values  for  the  trim 
settings  and  forces  and  moments  on  each  rotor 
during  the  tunnel  tests. 

EFFECT  OF  RESPONSE  ERROR 

The  blade  flexibility  is  handled  in  the  blade  loads  and  re- 
sponse program  through  the  use  of  normalized  rotating  blade 
mode  shapes  and  generalized  coordinates  which  determine  how 
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much  each  individual  mode  contributes  to  the  total  blade 
motions.  For  the  first  iteration,  no  blade  response  is  con- 
sidered (i.e.,  the  generalized  coordinates  are  all 
identically  zero  around  the  azimuth)  and  the  blade  loading 
which  would  occur  for  these  "rigid"  blades  is  computed.  The 
response  to  this  force  and  moment  distribution  is  then  de- 
termined by  computing  the  azimuthal  variation  which  each 
generalized  coordinate  (i.e.,  mode)  must  have  in  order  to 
make  the  blade  responses  compatible  with  these  forces.  The 
blade  response  has  a direct  effect  on  angle  of  attack  (and 
thus  blade  loading)  through  elastic  twist  and  an  indirect 
effect  through  plunging  velocity,  etc  . Thus,  the  blade 
loading  distribution  must  be  recomputed  with  the  blade  re- 
sponses taken  into  account.  Computing  the  blade  responses 
for  this  new  loading  then  provides  the  response  values  for 
the  second  iteration.  The  change  in  the  blade  responses 
between  successive  iterations  is  measured  by  the  response 
error,  e , which  is  defined  as 
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where  M is  the  number  of  modes  input  for  the  calculations, 
Na  is  the  number  of  azimuthal  steps  per  revolution,  and 

represents  the  generalized  coordinate  for  the  mode 
1 3 th 

at  the  i azimuth  station  with  the  superscript  2 designat- 
ing the  values  for  the  latest  iteration  and  a superscript 
of  1 for  the  values  on  the  previous  iteration.  When  the 
loading  and  the  responses  become  truly  compatible,  no 
changes  will  occur  in  the  loading  distribution  or  the 
blade  response  and  then  e will  drop  to  zero. 


A series  of  runs  was  made  for  th^  ABC  rotor,  with  smaller 
and  smaller  values  being  require*,  for  the  response  error. 
The  resulting  values  for  the  performance  parameters  and  for 
the  control  derivatives  are  shown  in  Figures  26  and  27. 

The  number  of  iterations  required  to  satisfy  the  various 
values  for  the  response  error  are  given  below: 
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response  error:  0.05  0.02  0,01  0.005 

number  of  iter- 
ations required:  7 9 10  13 

Good  values  were  obtained  for  all  of  the  performance  para- 
meters except  the  pitching  moment  with  er  = 0.05.  The  thrust 

and  yaw  control  derivatives  did  not  change  significantly  for 
0.02£er<^  0.05.  As  shown  in  Figure  27,  the  pitch  and  roll 

control  derivatives  did  not  vary  smoothly  with  er,  so  it  is 

unclear  what  response  error  is  required  to  provide  good  re- 
sults for  these  derivatives  in  this  particular  case.  Making 
er  = 0.02  would  provide  a fairly  good  representation  for  the 

pitch  and  roll  control  and  excellent  results  for  all  of  the 
other  variables  considered  here.  The  extra  iterations  re- 
quired to  make  the  response  error  less  than  0.02  are  not 
excessive  for  this  case,  so  a response  error  of  0.02  should 
be  used  in  future  runs  with  the  blade  loads  and  response 
program. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  main  conclusions  reached  from  this  study  are  that  wake 
effects  are  quite  important  to  the  performance  and  to  the 
stability  of  rotor  systems  which  (1)  have  overlapping  rotors 
so  that  the  wake  from  one  passes  through  the  other  or  (2) 
are  operating  at  low  speeds  so  that  the  wake  remains  nearby. 
More  detailed  conclusions  are: 

1.  Excellent  agreement  exists  between  the  measured 
and  calculated  thrusts  for  the  conventional  rotor 
when  the  measured  trim  settings  are  used. 

2.  A conventional  single  rotor  operating  at  moderate 
advance  ratios  and  shaft  tilt  angles  will  not  be 
sensitive  to  a deformed  wake.  However,  the  effects 
of  the  wake  are  significant  in  near  hover  condi- 
tions and  during  transition. 

3.  The  CTR  system  will  probably  always  be  sensitive 
to  the  v/ake.  The  torsional  moment  applied  by  the 
flap  is  affected  in  large  part  by  the  wake-induced 
velocity  distribution,  so  the  CTR  is  more  sensi- 
tive to  wake  effects  than  other  single-rotor 
systems  are. 

4.  The  wake  caused  significant  changes  in  the  per- 
formance and  stability  and  control  derivatives  on 
the  tandem-rotor  system  due  to  the  rotor  overlap. 

5.  The  counterrotating,  coaxial  rotor  system  of  the 
ABC  is  very  sensitive  to  wake  effects.  This  was 
especially  true  for  the  lower  rotor,  which  ap- 
proached or  went  into  autorotation  for  the  aft 
shaft-tilt  angles  considered  here  for  comparison 
with  experiment. 

6.  The  nonuniform  induced  velocity  distribution  on 
the  ABC  rotors  caused  a large  difference  to  exist 
between  the  results  with  uniform  inflow  and  those 
with  wake  effects  included  for  the  rolling  and 
pitching  moments.  The  net  rolling  and  pitching 
moments  usually  changed  direction  between  the  cases 
with  and  without  wake  effects  included. 

7.  A good  agreement  existed  between  the  trends  mea- 
sured for  the  values  of  the  stability  and  control 
derivatives  for  the  ABC  and  those  calculated  using 
wake-induced  velocities.  Better  agreement  might 
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be  obtained  for  the  cases  with  wake  effects 
included  if  the  rotor  were  retrimmed  with  the 
wake  included. 

Perturbing  the  wake  geometry,  and  consequently  the 
wake  influence  coefficients,  has  little  effect  on 
the  performance  parameters  for  the  ABC  in  forward 
flight.  The  performance  is  sensitive  to  wake  geo- 
metry variations  when  in  hover,  however. 

Most  of  the  variables  considered  herein  are  re- 
latively insensitive  to  blade  response  convergence. 
A response  error  of  about  0.02  produces  the  best 
results  overall  without  requiring  too  many  iter- 
ations. 


Based  on  the  conclusions  reached  in  this  investigation,  it 
is  recommended  that: 

1.  An  automatic  trim  procedure  be  set  up  using  the 
blade  loads  and  response  program,  especially  if 
rotors  are  to  be  trimmed  with  wake  effects  in- 
cluded. 

2.  A detailed  wake  and  blade  loading  investigation 
be  conducted  at  different  flight  conditions  for 
the  CTR  and  the  tandem-rotor  systems  in  order  to 
better  define  and  understand  the  wake  inter- 
actions and  its  effect  on  performance  and  stabil- 
ity. Of  particular  interest  would  be  the  reason 
for  the  power  increase  with  thrust  loss  on  the 
CTR  when  comparing  results  with  wake  effects  to 
those  with  uniform  inflow. 

3.  The  ABC  rotor  be  retrimmed  in  hover  and  in  forward 
flight  with  wake  effects  included.  The  stability 
and  control  derivatives  with  wake  effects  included 
should  then  be  computed  using  these  new  trim  set- 
tings for  their  basis. 

4.  The  actual  test  values  for  the  blade  pitch  settings 
on  each  rotor  and  the  individual  rotor  perfor- 
mance parameters  (if  such  are  available)  be  used 
for  a more  detailed  comparison  with  ABC  rotor  ex- 
perimental results. 

5.  In  considering  wake  effects,  the  wake  geometry 
and  wake  influence  coefficients  be  computed  for 
any  variations  in  flight  conditions  in  or  near 
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hover.  The  wake  geometry  for  a flight  condition 
for  y greater  than  about  0.2  is  probably  satis- 
factory for  other  similar  conditions. 

6.  A response  error  of  0.02  be  used  when  doing  the 
blade  loads  and  response  analysis. 

Further  development  should  be  conducted  in  order  to: 

1.  Include  the  full  aircraft  to  find  forces  and  mo- 
ments at  the  center  of  gravity. 

2.  Automatically  trim  an  entire  aircraft. 

3.  Calculate  dynamic  stability  derivatives  for  the 
full  aircraft. 


4.  Conduct  dynamic  stability  analyses  for  the  full 
aircraft. 
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The  angle  of  attack  for  zero  lift  varied  with  flap  deflection  for  the 
CTR  as  follows:  aT  = -1.0  -0.436 


Table  5.  Flight  Conditions  and  Trim  Settings 


upper  rotor  had  5 deg.  of  built  in  precone. 


Table  6.  Effect  of  Wake  on  Single  and  Tandem  Rotor  Performance  Parameters 
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Table  8.  Effect  of  Wake  on  Control  Derivatives  for  Single-Rotor  Systems 
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These  values  are  questionable  since  a small  "esponse  divergence  on  the 
Initial  Iterations  with  perturbed  pitch  settings  stopped  the  Iterations. 
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Effect  of  Wake  on  Control  Derivatives  for  Tandem  Rotor 
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Table  10.  Effect  of  Wake  on  ABC  Control  Derivatives  in  Hover 
for  0°  Shaft  Tilt  Angle 
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(2)  The  derivatives  with  respect  to  shaft  angle  were  not  computed  correctly  in  these  cases. 


•del  Showing  the  "Full  Mesh"  Wake  and 
idified"  Wake. 
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Figure  6.  Chordwise  Distribution  of  Bound  Circulation  and 
Downwash  for  an  Oscillating  Airfoil. 


Figure  7.  Transformation  of  Shears  and  Moments  on  Chordwise 
Axes  to  the  Rotor  Shaft  Axes. 
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Figure  16.  Effect  of  Wake  on  Thrust  and  Power  for  Each  Rotor 
of  the  ABC  in  Hover,  Ttotai  * 10,000  lb. 
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PRINCIPAL  DERIVATIVES 
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COUPLING  DERIVATIVES 

Figure  20.  Effect  of  Wake  on  the  Control  Derivatives  for 
a Tandem  Rotor,  u « 0.20. 
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Using  Wake  Geometry  Determined  For: 


Case  no. 
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Figure  22.  Effect  of  Different  Wake  Geometries  on  Control 
Derivatives  for  the  ABC  in  Hover,  Case  6, 
Reference  Thrust  = 15,000  Lb. 
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Figure  23.  Effect  of  Shaft  Angle  on  the  ABC  Stability  and 
Control  Derivatives,  p-0.208,  T.  . .=22,000  Lb, 


gj  Experimental  Results 
Q Uniform  Inflow 
YA  Wake- Induced  Inflow 


3ABi J i 5 


1 J 1 5 


3A6 oj  l 5 


CONTROL  DERIVATIVES 


LONGITUDINAL  STABILITY  DERIVATIVES 

The  value  for  each  derivative  when  T=22,000  lb  is  normalized 
using  its  value  when  T=1U,000  lb. 

Figure  24.  Effect  of  Thrust  on  the  ABC  stability  and  Control 
Derivatives.  u=0.208,  a_=4°. 


2200 


81 


Figure 


i w n wi"i( Tt. «. , iv m': nww«jv^.r u Wm WISPQIJ ■PPPW' ■ 1 ^ % W'SIM^OT  V P^tRlW  WWW.A'kUJ 


APPENDIX 

USERS  MANUAL  FOR  THE  WAKE  GEOMETRY  AND  BLADE  LOADS 
AND  RESPONSE  PROGRAMS  FOR  ADVANCED  ROTOR  SYSTEMS 

INTRODUCTION 


The  programs  to  be  described  in  this  users  manual  were  devel- 
oped in  order  to  determine  the  effect  of  the  rotor  wake  on 
the  performance  and  stability  characteristics  of  advanced 
rotor  systems.  These  rotor  systems  include  the  conventional 
type,  the  Telescoping  Rotor  Aircraft  system,  the  Controllable 
Twist  Rotor,  and  tandem  and  coaxial  counter-rotatTng  rotor 
systems.  For  dual-rotor  systems,  the  rotational  speed,  number 
of  blades  per  rotor,  and  the  number  of  airload  points  must  be 
the  same  for  each  rotor  but  all  other  blade  characteristics 
can  be  different  between  the  rotors.  The  wakes  shed  from 
both  rotors  will  fully  interact  with  each  other  when  the  wake 
geometry  calculations  are  done.  The  blades  can  be  canti- 
levered at  the  root,  articulated,  or  teetering.  The  mode 
shapes  and  frequencies  must  be  supplied  to  the  Blade  Loads 
and  Response  program.  These  shapes  and  frequencies  can  be 
obtained,  for  example,  using  the  Blade  Frequency  program  des- 
cribed by  S.  Gene  Sadler  in  NASA  CR-112071.* 


The  Wake  Geometry  program  determines  the  location  of  the  vor- 
tex elements  in  the  wake  and  their  effect  on  the  induced 
velocity  distribution  in  the  rotor  plane.  The  wake  is  self- 
deforming and  if  two  rotors  are  involved  their  wakes  will 
interact  with  and  deform  each  other.  The  radial  and  azimuth- 
al variations  £ the  bound  circulation  produces  the  shed  and 
trailing  vortices  which  are  in  the  wake.  The  bound  circu- 
lation strength  is  calculated  using  the  lift  distribution  on 
the  blades  and  this  lift  is  computed  using  linear  aerody- 
namics. Elastic  twist  and  blade  flapping  (up  to  2/rev)  can 
be  input  for  the  program  but  no  other  blade  responses  are 
considered.  The  wake-induced  velocity  coefficients  are  com- 
puted using  the  Biot-Savart  law  which  relates  the  location 
and  strength  of  the  vortex  elements  in  the  wake  to  the 
velocity  which  is  induced  at  the  blades.  The  use  cf  these 
coefficients  by  another  program  allows  the  bound  circulation 
distribution  to  be  changed  and  then  the  induced  velocity  dis- 
tribution to  be  adjusted  accordingly. 

The  response  of  the  blades  effects  their  angle-of-attack  dis- 


J 

Sadler,  S.  Gene:  INFORMAL  FINAL  REPORT  ON  BLADE  FREQUENCY 

PROGRAM  FOR  NONUNIFORM  HELICOPTER  ROTORS,  WITH  AUTOMATED  FRE- 
QUENCY SEARCH,  Rochester  Applied  Science  Associates,  NASA 
CR-112071  or  RASA  Report  No.  72-01,  April  1972. 
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tribution  and  thus  effects  the  loads  generated  by  the  blades 
which  in  turn  change  the  responses  of  the  blade.  Thus, 
iterations  are  carried  out  in  the  Blade  Loads  and  Response 
program  until  the  loads  and  response  are  compatible.  If  wake- 
induced  velocity  distributions  are  used,  then  the  bound 
circulations  are  calculated  from  the  blade  loading.  Multi- 
plying these  circulations  by  the  wake-induced  velocity 
coefficients  then  provides  the  induced  velocity  distribution 
which  can  be  used  in  recomputing  the  loads.  The  forces  and 
moments  transferred  to  the  hub(s)  by  all  of  the  blades  on 
each  rotor  are  calculated  and  these  hub  forces  and  moments 
are  then  transferred  to  the  C.G.  of  the  aircraft.  Perturbing 
the  shaft  angle  or  blade  pitch  settings  will  change  the 
forces  and  moments  at  the  C.G.  and  the  stability  and  control 
derivatives  are  computed  accordingly. 

The  hierarchy  charts  for  the  two  programs  are  shown  in  Fig- 
ures 28  and  29.  Descriptions  of  the  input  and  output  for  the 
two  programs  are  included  in  this  manual.  Sample  cases  for 
the  programs  are  provided  under  a separate  cover,  however, 
the  input  and  primary  output  for  one  such  case  is  shown  in 
Figures  30-34.  Tape  or  disk  unit  no.  4 must  be  set  up  to 
store  the  wake-induced  velocity  coefficients  from  the  Wake 
Geometry  program;  unit  no.  8 is  used  for  the  bound  circula- 
tions computed  by  the  Wake  Geometry  program  and  unit  no.  10 
is  used  to  store  the  variables  needed  for  restarting  the 
wake  geometry  calculations  if  an  earlier  run  does  not  finish. 
Units  4 and  8 are  used  by  the  Blade  Loads  and  Response  pro- 
gram only  if  wake-induced  velocity  distributions  are  to  be 
used  for  its  calculations. 
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Most  of  the  input  variables  for  the  wake  geometry  and  blade 
loads  and  response  programs  have  been  described  in  sufficient 
detail  in  the  preceding  tables.  An  extended  explanation  for 
some  of  those  variables  is  given  below. 

Linear  interpolation  is  used  to  compute  the  lift  curve  slope, 
zero  lift  angle,  ehordlength,  twist,  etc.,  from  the  distribu- 
tions input  to  the  wake  geomerry  program.  If  a distribution 
is  constant  or  linear  then  only  two  values  for  that  variable 
need  to  be  input:  the  value  at  the  root  and  the  value  at  the 
tip.  Nonlinear  distributions  must  be  broken  down  into  sections 
which  are  essentially  linear.  The  resulting  radial  positions 
and  the  value  of  the  variable  (e.g.,  twist)  at  those  radial 
positions  are  then  input  to  che  wake  geometry  program. 

The  number  of  blades  on  each  rotor  is  input  as  NIB  on  Card  #5 
for  the  wake  geometry  program  and  as  NB  on  Card  #5  for  the 
blade  loads  and  response  program.  The  programs  are  currently 
dimensioned  for  two  (2)  rotors  with  four  (4)  blades  each.  A 
single  rotor  with  six  (6)  or  eight  (8)  blades  can  be  modeled 
by  two  rotors  with  three  (3)  or  four  (4)  blades  each.  Both 
rotors  must  rotate  in  the  same  direction  (e.g.,  DIR  = 1.0  for 
both  rotors)  and,  with  PSIR  = 0.0  on  the  first  rotor,  PSIR 
will  be  60°  or  45°  for  the  second  rotor.  A rotor  with  five  or 
seven  blades  cannot  be  considered  since  the  spacings  between 
the  blades  will  not  be  equal.  Variable  geometry  rotors  can  be 
handled  this  same  way.  If,  for  example,  a four-bladed  rotor 
with  30°  blade  phasing  is  to  be  considered,  then  two  2- 
bladed  rotors  are  required  with  DIR  = 1.0  for  both  and  PSIR  » 
30.0  for  one  of  them. 

The  number  of  wake  points  to  be  considered  is  input  as  WW  on 
Card  #5  for  the  wake  geometry  program  and  as  NW  on  Card  #5  for 
the  blade  loads  program.  The  value  used  is  normally  WW  = 
NREV*NA  + 1 where  NREV  is  the  number  of  revolutions  saved  and 
NA  is  the  number  of  azimuthal  steps  per  revolution.  A frac- 
tional number  of  revolutions  can  be  used  if  some  other  value 
is  used  for  WW.  For  example,  if  NA  = 12  and  WW  * 17,  then 
1-1/3  revolutions  of  wake  are  used  to  compute  the  wake-induced 
velocities  at  the  blades.  The  extra  1 is  needed  since  the 
blade  itself  counts  as  wake  position  number  one. 

The  variables  NWKRQ  and  NUWKPT  on  Card  #7  of  the  wake  geometry 
program  are  nonzero  if  the  user  wants  to  compute  wake-induced 
velocities  at  arbitrary  points  in  space.  NUWKPT  is  the  number 
of  points  at  which  this  velocity  is  desired.  The  (x,y,z) 
coordinates  for  each  of  these  points  is  then  entered  on 
Card  #24. 

Blade  flapping  angles  (up  to  the  second  harmonic)  are  input  on 
Card  #20  of  the  wake  geometry  program  if  NCALB  is  zero  on 
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Card  #7.  If  NCALtt  = 1,  then  the  flapping  angles  are  computed 
internally  using  the  blade  mass,  inertia,  thrust,  etc.,  pro- 
vided on  Card  #21a-21h. 

The  reference  azimuth  position  is  input  as  PSIR  on  Card  #9  of 
the  wake  geometry  program  and  Card  #25  of  the  blade  loads  pro- 
gram. It  is  used  to  take  care  of  phasing  between  the  first 
and  second  rotors  and/or  a nonzero  azimuthal  starting  position 
for  either  rotor.  For  example,  a tandem  rotor  system  with 
four  blades  on  each  rotor  might  have  PSIR  = 0.0  for  the  forward 
rotor  and  PSIR  = 45.0°  for  the  aft  rotor.  PSIR  is  zero  when 
the  reference  blade  on  the  rotor  lies  downstream  and  is  in-line 
with  the  free  stream  velocity. 

The  lift,  drag,  and  moment  coefficients  can  be  obtained  in 
three  different  manners  in  the  blade  loads  program  depending  on 
the  value  input  for  I AERO.  If  I AERO  = 1 then  the  coefficients 
are  computed  internally  using  polynomials  which  have  been 
curve-fit  to  the  measured  data  for  the  NACA  0012  or  NACA  0015 
airfoils.  If  NAIR  = 0012  then  the  data  for  the  NACA  0012 
airfoil  is  used. 


If  I AERO  * 2 then  tables  for  the  incompressible  values  for  C^, 
C^,  and  C^  are  to  be  read  in  for  one  or  more  symmetric  airfoils. 

The  two-dimensional  incompressible  values  for  C^,  C^,  and  Cm  are 

input  on  one  card  for  each  angle  of  attack  (o  * 0°,  2°,  4°,  6°, 

. . . ,2* (IA-1) °) . This  block  of  data  is  then  repeated  with  the 
values  corresponding  to  each  successive  airfoil  type.  Suppose, 
for  example,  that  the  NACA  0006,  NACA  0009,  and  NACA  0012  air- 
foils are  to  be  used  with  the  angle  of  attack  ranging  up  to  16°. 
In  this  case,  IA  = 9 and  IT  = 3,  and  the  data  cards  would  appear 
as  follows: 


cl 

Cd 

cm 

(deg) 

Airfoil 

type 

0.00 

0.006 

0.00 

0 

0006 

0.21 

0.007 

0.00 

2 

0006 

0.41 

0.008 

0.00 

4 

0006 

0.61 

• 

0.015 

• 

0.00 

• 

6 

• 

0006 

• 

0.85 

0.20 

-.125 

1 1 

0006 

0.00 

0.0075 

0.00 

0 

0009 

0.20 

• 

0.0080 

• 

0.00 

• 

2 

• 

0009 

• 

1.06 

0.2 

-.080 

d 

0009 

0.00 

• 

0.008 

• 

0.00 

• 

0 

• 

0012 

• 

1.50 

.035 

-.005 

d 

0012 
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Note:  The  airfoil  type  and  o are  shown  on  the  previous  page 
only  for  the  sake  of  clarity.  They  may  be  put  on  the  cards 
if  the  user  so  desires,  but  they  are  not  required. 

Data  for  any  symmetric  airfoil  can  be  input  and  used  in  the 
blade  loads  program  as  it  now  stands.  To  make  the  values  input 
for  HAIR  in  this  case  correspond  to  the  code  numbers  for  such 
airfoils,  a few  changes  are  required  in  subroutine  CCC1. 

Suppose,  for  example,  that  the  symmetric  airfoil  types  to  be 
used  on  a given  rotor  blade  were  to  have  code  numbers  of 
65012,  65015,  and  65018.  Then  the  following  statements  at  the 
beginning  of  subroutine  CCCl, 

IS  = (ISP  - 3)/3 

IF  (IS  .GT.  ISMAX)  GO  TO  30950 

should  be  eliminated  and  replaced  by  the  following: 

IS  * 0 

IF  (ISP  .EQ. 65012)  IS  = 1 
IF  (ISP  .EQ. 65015)  IS  « 2 
IF  (ISP  .EQ. 65018)  IS  = 3 
IF  (IS  .EQ.  0 ) GO  TO  30950 

This  can  be  expanded  to  include  additional  airfoil  code  numbers 
(up  to  a maximum  of  7)  and  any  arbitrary  code  numbers  appear 
within  these  statements.  In  the  preceding  example,  the  values 
for  the  lift,  drag,  and  moment  coefficients  must  be  read  in  for 
a = Q® , 2®,  4®,  6®,  ...,2*(IA-1)®  for  airfoil  65012.  The  same 
tables  then  follow  with  the  values  for  airfoils  65015  and  65018. 

If  IAERO  = 3,  then  values  for  the  coefficients  at  various  angles 
of  attack,  Mach  number,  and  flap  deflection  angle  must  set  up 
for  input  with  the  tables  being  filled  in  that  order.  Nine  (9) 
values  are  read  on  each  card  so  that  the  lift,  drag,  and  moment 
coefficients  can  be  input  for  three  successive  angles  of  attack 
(at  the  current  Mach  number  and  flap  deflection  angle)  on  one 
card.  As  one  very  brief  example,  suppose  that  IA  = 5, 

IMOCK  = 2,  and  IFLP  = 3.  Then  the  data  cards  would  be  set  up 
as  shown  on  the  next  page. 
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WAKE  GEOMETRY  OUTPUT 


The  output  from  the  Wake  Geometry  program  can  be  broken  down 
into  four  sections  which  cover: 

1.  printout  of  the  input  data. 

2.  interpolated  values  for  the  blade  characteristics 
at  the  airload  points. 

3.  intermediate  values  for  the  blade  loading  and 
circulations  as  the  wake  is  developed. 

4.  final  values  for  the  wake-induced  velocity  in- 
fluence coefficients  and  the  wake  geometry. 

INPUT  DATA 

The  input  data  is  presented  in  the  same  order  as  it  appeared 
on  the  cards  for  input  to  the  computer.  The  one  exception  to 
this  is  the  values  for  the  RCAPS  if  they  are  read  in  rather 
than  computed  internally.  The  description  for  the  input  data 
is  given  in  the  preceding  section  of  this  documentation.  An 
example  of  the  data  for  the  ABC  rotor  is  shown  in  Figure  32. 

INTERPOLATED  VALUES 

RCAP  is  the  radial  coordinate  of  the  trailing  vortices  given 
in  units  of  "feet".  These  values  are  computed  internally  if 
ISWRAD-0  and  were  read  in  if  ISWRAD=1. 

RBL  represents  the  radial  locations  of  the  airload  points. 

These  values  are  nondimensional  and  are  defined  by 

RBLi  * (RCAPi  + RCAPi+1)/2R 
where  R is  the  rotor  radius. 

The  blade  characteristics  at  each  RBL  are  obtained  by  inter- 
polation from  the  input  data  at  its  arbitrary  radial  locations. 
The  characteristics  involved  are: 

chordlength,  ft 

lift  curve  slope,  dc^/da,  per  radian 

zero  lift  angle  with  no  flap  deflection,  (0) , deg 

o 

change  in  zero-lift  angle  with  flap  deflection, 

daLo/dfif , deg/deg 

built-in  twist,  <Mr) , positive  nose  up,  deg 
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An  example  of  the  interpolated  values  for  the  ABC  rotor  is 
shown  in  Figure  33a. 

INTERMEDIATE  CALCULATIONS 

One  set  of  values  during  the  intermediate  calculations  for 
the  ABC  rotor  is  given  in  Figure  33b. 

For  each  azimuthal  step  taken  during  the  build-up  of  the  shed 
wake,  the  following  data  is  printed  out: 


r(l> 

i 

for  i=l,2,..., 

NTV1*NR0T*NIB 

(2) 
r . 

X 

for  x—1 ,2 . . » . , 

NTV1*NR0T*NIB 

ITR  GTE ST 

MSET 

( 3\ 

XK  and  r}  ’ 

for  i— 1 , 2 , . . . , 

NTVl*NROT*NIB 

L‘i 

for  i=l, 2 , , 

NTV1*NR0T*NIB 

NAS  NW 

NWSTRE 

These  variables  are  defined  as  follows: 

r±  is  the  nondimensional  bound  circulation  at  each  airload 

point  on  each  blade  of  each  rotor  for  the  current  azimuthal 
position.  This  circulation  is  calculated  from  the  relation 

Ti  * 1/2  c C j a /u2+V2  /(fiR2) 
a 


: with  the  angle  of  attack, a,  calculated  under  three  different 

conditions.  These  conditions  correspond  to  the  superscripts 
( for  r.  and  involve  various  assumptions  on  the  induced  velo- 

> A 

) city  distribution.  These  assumptions  are: 

s 1)  is  identically  zero  at  all  points. 

? 2)  v^  is  computed  from  the  wake-induced  velocity  in- 

: fluence  coefficients  and  bound  circulation  values 

computed  for  the  azimuthal  steps  up  to  and  includ- 
{ ing  the  last  one. 

3)  v^  is  computed  as  in  step  2)  except  that  the  current 

k azimuthal  step  is  also  included. 
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In  step  3)  above,  ri  changes  the  induced-velocity  dis- 
tribution and  thus  the  angle  of  attack  by  which  it  is 
computed.  Thus,  iterations  must  be  carried  out  to  make  the 
bound  circulations  and  the  wake- induced  velocities  correspond 
to  each  other.  These  iterations  involve  ITR,  GTEST,  MSET, 
and  XK  where 

ITR  is  the  number  of  iterations  required  for  con- 
vergence of  the  r^s. 

GTEST  is  a measure  of  the  differences  for  the  rs 
between  the  last  two  iterations. 

MSET  is  the  blade  position  index  (see  Figure  35) . 

XK  is  a measure  of  the  induced  velocity  at  the 
last  airload  point  on  the  last  blade.  It  should 
be  ignored. 

L^  is  the  lift  per  unit  length  at  each  airload 

point  on  each  blade  of  each  rotor  with  units  of 
lb/ft. 

NAS  is  the  total  number  of  azimuthal  steps  taken. 

NW  is  the  number  of  steps  of  wake  in  the  full 
mesh  and  will  invariably  equal  NANRM. 

NWSTRE  is  the  number  of  steps  included  in  the  full 
and  modified  mesh  (with  the  blade  itself  counting 
as  1) . 

NWSTRE=NAS  until  the  wake  is  fully  developed  and 
computations  begin  for  the  final  values. 

NT VI  is  the  number  of  airload  points  (i.e.,  number 
of  trailing  vortices  minus  one) . 

NROT  is  the  number  of  rotors. 

NIB  is  the  number  of  blades  per  rotor. 

FINAL  VALUES 

A portion  of  the  final  values  for  the  ABC  rotor  is  shown  in 
Figure  33 c.  a brief  description  of  these  blocks  of  numbers 
follows. 
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The  wake-induced  velocity  influence  coefficients  and  the  posi- 
tions and  velocities  of  the  end  points  of  each  vortex  element 
in  the  wake  are  given  in  the  final  output.  The  number  of 
blocks  of  this  data  which  appear  is  equal  to  NA/NIB.  Data 
is  presented  for  each  blade  of  each  rotor  (i.e.,  NIB  sets) 
within  each  block  and  so  the  full  azimuth  is  covered. 

The  first  set  of  data  to  appear  in  each  block  for  the  final 
values  is  a set  of  tables  for  the  wake-induced  velocity  co- 
efficients for  each  airload  point  on  each  blade  of  each  rotor 
at  the  azimuthal  location  corresponding  to  the  total  block. 
Each  table  is  headed  by  the  integer  MSET  which  identifies  the 
airload  point  which  corresponds  to  the  given  table  of  coeffi- 
cients. The  indexing  scheme  is  illustrated  in  Figure  35. 

The  number  of  entries  in  each  table  is  NTV1*NR0T*NA  which  is 
equal  to  the  total  number  of  bound  circulations  which  would 
influence  the  downwash  at  the  airload  point.  These  numbers 
are  printed  out  in  the  same  sequence  in  which  they  are  stored 
on  disk  or  tape  for  use  by  the  Blade  Loads  and  Response  pro- 
gram. 

The  intermediate  calculations  described  earlier  are  also  done 
for  the  final  values.  The  bound  circulations  and  the  lift 
distributions  are  presented  following  the  wake-induced  veloc- 
ity coefficients  just  as  they  were  for  the  intermediate 
calculations. 

The  next  set  of  data  to  appear  is  the  wake-induced  velocity, 
nondimensionalized  by  the  rotor  tip  speed,  in  each  direction 
at  the  end  points  of  each  vortex  element  in  the  full  mesh. 
These  velocity  arrays  are  VX(i,j),  VY(i,j),  and  VZ(i,j)  with 
"iM  varying  from  1 to  NA  and  " j"  varying  from  1 to  NTV*NIB* 
NROT.  The  subscript  "i"  refers  to  azimuthal  positions  in  the 
full  mesh  with  the  blade  itself  counting  as  1 . The  subscript 
"j"  covers  each  trailing  vortex  on  each  blade  of  each  rotor. 

The  spatial  location  of  the  end  point  of  each  vortex  element 
in  the  full  mesh  follows  the  velocity  output.  The  indices 
i and  j cover  the  same  range  and  have  uhe  same  meaning  that 
they  did  for  the  velocity  arrays. 

The  induced  velocity  and  the  position  for  the  end  points  of 
the  vortex  elements  in  the  modified  mesh  conclude  the  data 
given  in  each  block  of  the  final  values.  The  index  j varies 
from  1 to  NTVM*NIB*NuOT  and  so  covers  the  number  of  trailing 
vortices  in  the  modified  mesh  for  each  blade  of  each  rotor. 
The  index  i varies  from  1 to  NA*NREV  and  so  covers  each  azi- 
muthal step  for  the  number  of  revolutions  (NREV)  saved  in 
the  wake.  The  first  NANRM-1  rows  represent  points  in  the 
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£ull  mesh  which  are  irrelevant  and  therefore  not  considered 
in  the  calculations.  Any  numbers  in  those  rows  are  to  be  ig- 
nored by  the  user.  It  might  be  noted  t*u  t the  location  of 
the  tip  vortex  is  the  same  for  the  full  wake  and  the  modified 
wake  at  the  boundary  between  the  two  meshes.  This  boundary 
exists  for  i=NANRM,  and  the  tip  vortex  is  located  at  a) 

J=n*NTV  for  the  full  mesh  and  b)  j=n*NTVM  for  the  modified 
mesh  with  n=l,2, . . ..  ,NIB*NROT. 

In  the  final  block  a full  table  of  the  bound  circulations  is 
printed  out  preceding  the  velocity  matrices.  The  values  in 
this  table  are  in  the  same  sequence  with  which  they  are  read 
onto  disk  or  tape  for  use  by  the  Blade  Loads  and  Response 
program  a a initial  estimates  for  the  distribution  of  the 
bound  circulations  with  no  blade  response. 


BLADE  LOADS  AND  RESPONSE  OUTPUT 

The  four  main  sections  of  the  output  from  the  Blade  Loads  and 
Response  program  are  concerned  with: 

1 . the  input  data . 

2.  the  intermediate  values  of  loading  and  blade 
response  which  are  computed  during  the  iterations 
required  to  make  the  loads  and  response  be 
compatible. 

3.  the  aerodynamic  loading  and  the  blade  responses 
for  each  rotor  after  the  iterations  converge  as 
well  as  the  forces  and  moments  transferred  to  the 
hub  of  each  rotor. 

4.  the  total  forces  and  moments  at  the  C.G.  and  the 
derivatives  with  respect  to  shaft  angle  or  blade 
pitch  settings  after  they  have  been  perturbed. 

INPUT  DATA 

The  general  input  data  appears  just  as  it  does  for  input 
with  the  cards  and  the  description  for  each  variable  is 
given  in  the  section  for  input  data.  The  input  value  for  IN, 
which  controls  the  unit  on  which  the  aerodynamic  and  mode 
shape  data  is  read  in,  is  omitted  in  this  printout,  however. 

A sample  set  of  the  input  data  for  the  ABC  rotor  is  given  m 
Figure  34a. 

The  tables  for  the  aerodynamic  coefficients  are  printed  out 
if  I AERO  is  2 or  3 and  if  NPRNT=2. 

The  rotor  properties  and  the  blade  properties  are  printed 
out  just  as  they  are  supplied  for  i:.put. 

11? 


The  torsional  damping  matrix  is  computed  from  input  values  as 
follows: 

SIGKJ  (K,J)  = DAMPC*  APHI (1,K) *APHI (1, J)  for  K*J 

and 

SIGKJ  (K,K)  = -2*SIG(K) *OMEGA(K) 

where 

DAMPC= torsional  damping  coefficient,  ftrlb-sec 
APHI (1 ,K)  = torsional  response  of  the  Kz  mode  at  the 
outboard  end  of  the  first  blade  element,  rad 

SIG^K)  = input  value  for  the  modal  damping  coefficient 
of  the  K^mode 
and 

OMEGA (K)  = natural  frequency  of  the  a mode,  rad/sec 

The  data  for  the  normal  modes  is  provided  by  input. 

The  values  are  those  supplied  as  input  except  that  the  sign 
of  w and  are  changed.  The  variables  defining  the  shape  of 

each  mode  are  v,w,  <|>,SI,  and  THETA  where 

v = flatwise  deflection,  positive  down,  ft 
w = edgewise  deflection,  positive  for  lead,  ft 
= torsional  deflection,  positive  nose  up,  rad 
SI  = 3v/3r 

and 

THETA  = 3w/3r 

The  frequency  for  each  mode  is  self-evident.  The  damp- 
ing term  is 

DAMPING  SIGMA  = SIG(K)  + DAMPC* [APHI(1,K) 3 2 


The  shear  and  moment  distribution  for  each  mode  concludes  the 
input  data.  The  flatwise  and  edgewise  shears  and  moments 
(MZ,-VY,MY,  and  VZ)  are  illustrated  in  Figure  7.  The  tor- 
sional moment,  T,  is  positive  nose  up. 

INTERMEDIATE  VALUES 


The  amount  of  information  which  is  printed  out  during  the 
iterations  varies  with  the  control  variable  NPRNT.  Very 
little  printout  is  produced  if  NPRNT=-2.  The  output  to  be 
described  in  this  section  is  that  which  occurs  for  NPRNT=1 
plus  the  more  significant  terms  among  those  additional  ones 
which  would  appear  for  NPRNT=2. 

The  local  airspeed,  angle  of  attack,  wake-induced  downwash, 
and  bound  circulation  are  given  for  each  airload  point  on 
each  rotor  at  each  azimuthal  step.  The  angle  of  attack  is 
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positive  nose  up,  the  downwash  is  positive  down,  and  the 
bound  circulation  will  be  positive  if  the  lift  is  positive 
(i.e.,  up).  There  will  be  NR1*NR0T  columns  in  each  table  to 
cover  the  airload  points  (as  read  in  through  the  RBL's)  for 
each  rotor.  The  first  NR1  columns  are  for  the  airload 
points  (RBL's)  on  the  first  rotor  and,  if  a second  rotor 
exists,  the  next  NR1  columns  will  be  for  the  second  rotor. 

The  azimuthal  position  of  the  blade  on  each  rotor  varies  from 
PSIR  to  PSIR+DIR* (NA-1) *360°/NA  for  the  rows  in  each  table. 

The  reference  azimuth  is  given  on  input  by  PSIR,  the  direc- 
tion of  rotation  is  counterclockwise  for  DIR=1.  and  clockwise 
for  DIR=-1,  and  the  number  of  azimuthal  steps  is  NA*  PSIR 
and  DIR  can  be  different  for  each  rotor  if  two  rotors  are  in- 
volved. This  tabular  scheme  for  presenting  the  radial  and 
azimuthal  variations  of  the  output  variables  for  each  rotor 
is  used  whenever  values  at  the  airload  points  are  to  be 
printed. 

IT1,IT2,IT3,  and  the  ERRORS  assocated  with  them  are  printed 
out  so  that  the  iterations  during  the  calculations  and  their 
rate  of  convergence  can  be  monitored.  The  innermost  itera- 
tion is  an  iterative  scheme  used  to  solve  a set  of 
simultaneous  equations  for  the  bound  circulations.  The  wake- 
induced  velocity  is  not  allowed  to  vary  within  this  iteration. 
The  iterations  on  the  innermost  loop  will  continue  until 
ERROR  <ALL1  or  IT1=NIT1.  For  the  second  loop,  the  wake-in- 
duced velocity  is  recomputed  to  correspond  to  the  new  set  of 
bound  circulations  and  bound  circulations  are  recomputed  to 
correspond  to  the  updated  velocity  distribution.  The  changes 
in  the  circulations  due  to  induced-velocity  corrections  then 
make  up  the  error  associated  with  the  second  loop.  Itera- 
tions will  continue  on  the  second  level  until  the  errct  on 
this  level  is  less  than  ALLl  or  until  IT2=NIT2.  Finally, 
the  revised  bound  circulations  represent  new  loading  distri- 
butions for  the  blades  and  so  the  blade  responses  must  be 
corrected  to  make  them  correlate  with  the  new  loading.  The 
changes  in  the  response  quantities  make  up  the  RESPONSE  ERROR 
which  is  associated  with  IT3.  Response  iterations  will  con- 
tinue until  RESPONSE  ERROR  <ALL'2  or  until  IT3  * NIT3. 

The  azimuthal  variation  of  the  generalized  force  for  each 
mode  is  printed  out  on  the  lines  labeled  FORC.  These  forces 
have  been  integrated  over  the  radius  of  the  blade  for  each 
azimuthal  position  from  PSIR  to  DIR*3<'.0°  (NA-1) /NA+PSIR. 

If  NPRNT=2 , a considerable  amount  of  additional  data  is 
printed  out  during  the  iterations.  The  more  significant 
portions  of  this  printout  include  1)  the  aerodynamic  forces 
acting  at  the  airload  points,  2)  the  shears  and  moments  at 
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the  mass  points  due  to  the  aerodynamic  loading,  3)  the  gener- 
alized coordinates,  and  4)  the  plunging  velocity  and  elastic 
twist  of  the  blade  at  each  airload  point. 

The  table  headed  LOADS  ON  BLADE  contains  the  azimuthal  varia- 
tion of  the  aerodynamic  loading  at  each  air lead  point.  FORCE 
Z is  in  the  direction  of  the  shaft  and  is  positive  up.  FORCE 
X is  in  the  plane  normal  to  the  shaft  and  is  positive  back 
(i.e.,  causing  blade  lag).  Both  of  these  forces  are  nondi- 
mensionalized  by  the  factor  pfi2R2(c/2).  The  moment  is 
measured  about  the  midchord,  is  positive  nose  up,  and  is 
nondimensionalized  by  the  factor  pft2R2(c/2)2.  The  printout 
for  the  last  iteration  is  dimensionalized  with  FORCE  Z and 
FORCE  X having  units  of  lb/ft  and  MOMENT  having  the  units  of 
ft-lb/ft. 

The  aerodynamic  loading  is  integrated  between  each  mass 
point  assuming  a linear  variation  between  airload  points,  and 
the  resulting  shears  and  moment  are  given  in  the  tables 
labeled  FV,  FW,  and  EMOME.  The  values  are  given  for  each 
mass  point  from  the  one  nearest  the  root  to  the  one  nearest 
the  tip.  FV  is  the  shear  normal  to  the  chord,  is  positive 
down,  and  is  given  with  units  of  lb.  FW  is  the  shear  along 
the  chordline  in  lb  and  is  positive  toward  the  trailing  edge. 
EMOME  is  the  torsional  moment  about  the  elastic  axis  at  the 
mass  points,  is  positive  nose  up,  and  has  units  of  ft-lb. 

The  generalized  coordinate  for  each  mode  and  its  first  two 
time  derivates  are  given  in  the  tables  labeled  CSI,  CSIDT, 
and  CS2DT.  The  columns  correspond  to  the  various  modes  and 
] the  rows  cover  the  azimuth  positions  from  PSIR  to  DIR*360° 

I + PSIR. 

The  plunging  velocity  and  elastic  twist  of  the  blades  at  the 
S airload  points  are  given  in  tables  labeled  by  VERTICAL 

? VELOCITY  AT  AIRLOAD  POINTS,  ROOT  PITCH  ANGLE  PLUS  RADIAL 

| ELASTIC  TWIST,  and  HARMONIC  ANALYSIS  OF  ELASTIC  TWIST.  These 

* tables  also  appear  in  the  final  output  for  each  rotor  and 

will  be  described  in  the  next  section. 

FINAL  OUTPUT  FOR  EACH  ROTOR 

4 

The  azimuthal  variation  of  the  aerodynamic  forces  and  tor- 
sional moment  at  each  airload  point  for  each  rotor  are 
printed  out  as  the  first  part  of  the  final  output.  The 
first  NRl  columns  correspond  to  the  airload  points  of  the 
first  rotor  and  the  right-hand  block  is  for  the  second  rotor, 
if  it  exists.  The  azimuth  varies  from  PSIR  to  PSIR  + DIR* 
360°  (NA-D/NA  for  each  rotor.  The  force  in  the  z-direction 
is  in  the  direction  of  the  shaft  and  is  positive  up. 
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The  force  in  the  x-direction  is  in  the  plane  normal  to  the 
shaft  and  is  positive  toward  the  trailing  edge  of  the  blade. 
These  forces  are  given  with  units  of  lb/ft.  The  torsional 
moment  is  referenced  to  the  aerodynamic  center,  is  positive 
nose  up,  and  has  units  of  ft-lb/ft. 

Each  rotor  is  treated  individually  following  the  output  for 
the  aerodynamic  forces  and  moments.  Blade  response  velocities 
and  slopes  are  given  for  each  mass  point  along  the  blade.  The 
flatwise  velocity  is  normal  to  the  chord,  positive  down,  and 
has  units  of  ft/sec.  The  chordwise  velocity  is  positive  to- 
ward the  leading  edge  in  the  direction  of  the  chord  and  has 
units  of  ft/sec.  The  torsional  deflection  angular  rate  is 
positive  nose  up  with  units  of  rad/sec. 

The  table  labeled  VERTICAL  VELOCITY  AT  AIRLOAD  POINTS  shows 
the  blade  motions  in  the  direction  of  the  shaft  due  to  the 
blade  response.  This  velocity  is  positive  up  and  is  nondi- 
mensionalized  by  the  rotor  tip  speed. 

The  ROOT  PITCH  ANGLE  PLUS  RADIAL  ELASTIC  TWIST  is  given  in 
radians  at  the  airload  points  of  each  rotor  and  is  positive 
nose  up.  This  angle  represents  the  geometric  angle  of  attack 
with  the  built-in  twist  removed  and  is  equal  to 

eQ  - A^cosiji  - B^sin^  + 4>e(ri,^) 

A harmonic  analysis  of  the  elastic  twist  is  carried  out  and 
the  Fourier  coefficients,  in  degrees,  for  the  twist  at  each 
airload  point  are  printed.  The  Fourier  series  for  the 
elastic  twist  can  be  written  as 

(NA/2) -1 

4>e(r.,*)  = aQ(ri)  + l [a^r^  cos  n*  + W sin 

n=l 


with  r^  varying  from  the  airload  point  nearest  the  root  to 
the  airload  point  nearest  the  tip. 

The  displacements  of  the  blade  at  each  mass  point  are  given 
as  a function  of  azimuth.  The  torsional  deflection  angle  is 
due  to  the  elastic  twist  of  the  blade;  it  is  given  in 
radians  and  is  positive  nose  up.  The  flatwise  displacement 
is  positive  nose  up.  The  flatwise  displaement  is  positive 
down,  with  units  of  ft,  and  is  measured  normal  to  the  chord. 
The  chordwise  displacement  is  positive  toward  the  leading 
edge  of  the  blade,  has  units  of  ft,  and  is  measured  parallel 
to  the  chord. 
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The  shears  and  moments  at  the  inboard  side  of  each  mass  point 
are  shown  in  their  respective  tables  as  they  vary  with  azi- 
muth. The  torque  is  positive  nose  up  with  units  of  ft-lb. 

The  units  for  the  y-moment  and  z-moment  are  ft-lb, and  the 
units  for  the  y-shear  and  z-shear  are  lb.  The  sign  conven- 
tions and  reference  axes  for  these  shears  and  moments  are  shown 
in  Figure  7.  A harmonic  analysis  is  conducted  for  the 
shears  and  moments  at  the  root, and  the  Fourier  coefficients 
for  these  root  values  are  given  at  the  head  of  each  table. 

These  coefficients  are  used  to  determine  the  forces  and  mo- 
ments which  are  transferred  to  the  hub. 

The  hub  drag  is  the  force  in  the  shaft  plane,  lb,  and  is  posi- 
tive aft.  The  thrust  is  the  force  in  the  direction  of  the 
shaft,  positive  up,  with  units  of  lb.  The  torque  is  a measure 
of  the  power  required  to  turn  the  rotor  and  has  units  of  ft-lb. 
The  pitching  moment  at  the  hub  is  positive  if  it  causes  the 
helicopter  to  pitch  nose  up  and  is  given  with  units  of  ft-lb. 
The  rolling  moment  at  the  hub  is  positive  if  it  causes  a 
right  roll  for  the  helicopter  and  has  units  of  ft-lb. 

VALUES  AT  THE  C.G.  AND  DERIVATIVES 

The  thrust  is  the  total  vertical  force  in  lb,  and  positive  up. 
The  "drag"  is  the  horizontal  force  in  lb  due  to  the  rotors 
and  is  positive  in  the  direction  of  flight.  This  horizontal 
force  includes  the  horizontal  component  of  the  thrust  from 
each  rotor.  The  TORQUE  is  the  yawing  moment  about  the  C.G. 
in  ft-lb  and  is  positive  if  it  causes  the  helicopter  nose  to 
yaw  to  the  right.  This  yawing  moment  includes  the  horizon- 
tal force (s)  at  the  rotor  hubs  times  their  lateral  offsets 
from  the  C.G.  The  horsepower  is  the  power  required  to  drive 
the  rotors.  The  pitching  moment  in  ft-lb  is  positive  nose  up 
and  includes  the  vertical  and  horizontal  force (s)  at  the 
rotor  hubs  times  their  horizontal  and  vertical  offsets  from 
the  C.G.  The  rolling  moment  in  ft-lb  causes  a right  roll  if 
it  is  positive  and  includes  the  vertical  rotor  forces  times 
their  lateral  offsets  from  the  C.G. 

If  the  shaft  angle  or  blade  pitch  settings  have  been  per- 
turbed (i.e.,NIP>l)  then  the  stability  or  control  derivatives 
are  computed  and  printed  out  following  the  summary  of  forces 
and  moments  for  the  perturbed  settings  about  the  C.G.  The 
derivatives  for  any  forces  have  units  of  lb/deg,  the  deriva- 
tives of  the  moments  are  in  ft-lb/deg,  and  the  variation  of 
power  required  with  the  shaft  angle  or  blade  pitch  setting 
is  in  HP/deg.  The  derivatives  which  are  computed  appear  in 
the  following  order: 
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Derivative  taken  with  respect  to: 
shaft  angle 
collective  pitch 
longitudinal  cyclic  pitch, 

lateral  cyclic  pitch,  A^ 

differential  collective  pitch,  A6  (for  two 
rotors)  0 

differential  longitudinal  cyclic  pitch,  AB^ 

(for  two  rotors) 


A portion  of  the  intermediate  and  final  values  for  each 
rotor  and  the  values  at  the  C.G.  for  the  sample  case  for  the 
ABC  rotor  are  shown  in  Figure  34b.  The  derivatives  with 
respect  to  collective  pitch  are  given  in  Figure  34c. 
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Figure  31.  Sample  Portion  of  the  Blade  Loads  and  Response  Input  for  the  CTR  Rotor 
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Figure  32.  Wake  Geometry  Input  Data  for  a Sample  Case  for  the  ABC  Rotor  System 
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Figure  33.  Continued. 
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Figure  34.  Continued 
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LIST  OF  SYMBOLS 


lateral  cyclic  pitch,  deg 


q 

,(k) 


*T 


(K) 


A(K),  a(k), 

a(k)  a(k) 
A$  • % ' 

.(K) 


mode  shape  quantities  representing  the  "q"  type 
mode  variable 

mode  shape  quantity  for  the  "qn  type  of  elastic 
deformation,  where  q = v,w,^,¥,  or  ^ 

coefficient  of  torsional  moment  at  ith  radial 
station  in  Kth  mode,  ft-lb 

mode  shape  quantities  of  the  Kth  mode  represent- 
ing the  lir.v  r flap  deflection,  linear  lead- lag 
deflection,  angular  torsional  deflection,  angular 
lead-lag  deflection  and  flapwise  bending  slope, 
respectively.  Linear  deflection  in  ft,  angular 
deflection  in  rad 


speed  of  sound,  ft/sec 


ao 


Vbn 


B, 


coefficient  of  zeroth  harmonic  of  elastic  twist 
with  respect  to  ip,  rad 

harmonic  coefficients  of  elastic  twist  with 
respect  to  <!<,  rad 

longitudinal  cyclic  pitch,  deg 


bi 


c 

c- 


blade  semichord  at  radius  r.,  nondimensionalized 
by  R 1 

blade  chord  length,  ft 

drag  coefficient 


torsional  damping  coefficient 


lift  coefficient 


lift  curve  slope 
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*at»v 


v . i,.;.’.vi''Jvn,)’.v.l.7‘"»<. :.'.  •'■p-'i?.';- «.  ,,;Ti-v?"  z-  t* 


HP 

io 

l 

M 

M 

M 

o 

My 

\ 

Mz 

M_ 


Mz  'Mz 
C1  S1 


horizontal  distance  from  elastic  axis  to  pitch 
axis,  positive  for  pitch  axis  ahead  of  elastic 
axis,  ft 

plunging  velocity  for  blade  section,  positive 
down,  ft/sec 

horsepower 

blade  element  torsional  mass  moment  of  inertia 
about  the  elastic  axis,  lb-ft-sec^ 

lift  per  unit  span,  lb/ft 

Mach  number  (=local  aero,  velocity/speed  of  sound) 
number  of  modes 

aerodynamic  moment  about  the  midchord,  positive 
nose  up,  ft-lb 

in-plane  bending  moment,  normal  to  rotor  plane, 
positive  down,  ft-lb 

chordwise  bending  moment,  normal  to  chordline, 
positive  down,  ft-lb 

out-of-plane  bending  moment,  parallel  to  rotor 
plane,  positive  toward  leading  edge  (i.e.,  blade 
bends  down  with  increasing  radius),  ft-lb 

flatwise  bending  moment,  in-line  with  chord, 
positive  toward  leading  edge,  ft-lb 

first  harmonic  cosine  and  sine  components  of 
pitching  moment  of  all  blades,  ft-lb 


M, 


M,  'M<* 

*c  *s 
C1  S1 


lumped  aerodynamic  twisting  moment  about  elastic 
axis,  lb-ft 

first  harmonic  cosine  and  sine  components  of 
rolling  moment  of  all  blades,  ft-lb 


m 

mo 


blade  element  mass,  lb-sec  /ft 

aerodynamic  moment  applied  at  a blade  station 
about  the  quarter  chord,  ft-lb 
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"T5*- * "'j 


number  of  azimuthal  steps  per  revolution 
number  of  radial  load  points  per  blade 
number  of  blades  per  rotor 

pitching  moment  at  rotor  hub,  positive  nose  up,  ft- lb 
rotor  torque,  ft-lb 

forcing  functions  corresponding  to  linear  flap, 
linear  lead-lag,  angular  twist,  angular  lead-lag 
motions,  and  flapwise  bending  slope  used  in 
computing  generalized  forces 

vortex  induced  velocity  at  a point  located  by  the 
vector  s,  ft/sec 

induced  velocity  at  a point  in  space  due  to  a 
straight  vortex  element,  ft/sec 

rotor  radius , ft 

rolling  moment  at  rotor  hub,  positive  for  right 
roll,  lb 

radial  distance  from  rotor  hub,  ft 
lateral  lift  offset,  ft 


position  vector  of  vortex  element 

position  vector  of  point  at  which  induced 
velocity  is  computed 

position  vector, 
thrust,  lb 

torsional  moment,  ft-lb 


tail  rotor  thrust,  lb 


thrust  of  upper  and  lower  rotors,  respectively, 
lb 
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tangential  velocity#  ft/sec 


local  airspeed 


. = /u2  + V2)  , 


ft/sec 


velocity  normal  to  rotor  plane,  positive  up,  ft/sec 
flight  velocity  , ft/sec 

shear  normal  to  the  rotor  plane,  positive  up,  lb 

flatwise  shear,  normal  to  chord,  positive  up,  lb 

vertical  zeroth  harmonic  component  of  shear  of 
all  blades,  lb 


V ,V 
yc,  ySl 


vertical  first  cosine  and  sine  harmonic  components 
of  shear  of  all  blades,  lb 


in-plane  shear,  parallel  to  rotor  plane,  positive 
toward  leading  edge,  lb 

chordwise  shear,  parallel  to  chord,  positive 
toward  leading  edge,  lb 

in-plane  first  harmonic  sine  component  of  shear 
of  all  blades,  lb 

flatwise  deflection,  positive  down,  ft 

the  chordwise  average  downwash,  ft/sec 

downwash  induced  by  the  wake,  positive  down,  ft/sec 

edgewise  deflection,  positive  toward  trailing 
edge,  ft 

climb  velocity,  positive  up,  ft/sec 


X,Y,Z 


root 


rectangular  coordinates  as  defined  in  Figure  3 
offset  of  flapping  hinge  from  hub,  ft 

yawing  moment,  positive  for  nose  right,  ft-lb 

distance  of  elastic  axis  forward  of  quarter  chord, 
ft 
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J.’"T 


a 

a. 


aL  (0) 


8 


30rBl,3Z 


r 


y(r, ?,<!») 


angle  of  attack#  rad 

forward  tilt  of  the  rotor  plane  with  respect  to 
the  shaft  axis  due  to  flapping,  rad 

geometric  pitch  angle  for  the  blade  section,  rad 
angle  of  attack  for  zero  lift,  rad 

angle  of  attack  for  zero  lift  with  no  flap 
deflection,  rad 

shaft  tilt  angle,  positive  aft,  rad 

blade  flapping  angle,  (6  = B + Bc  sin\|>  + B cos*), 
rad  ° c 

steady  and  first  harmonic  lateral  and  longitudinal 
blade  flapping  components,  rad 

2 

vortex  element  circulation,  ft  /sec 

2 

bound  circulation  strength,  ft  /sec 

chordwise  elemental  circulation  strength 
incremental  change  in  yaw  moment,  ft-lb 

incremental  change  in  rolling  moment,  ft-lb 


At 


6 ( A 6 ) 
o 


0 


time  increment,  sec 

differential  collective  pitch  angle  between  two 
rotor  systems,  rad 

azimuthal  increment,  deg 

flap  deflection,  positive  down  (6^  = 6q  6ssint|» 

+ 6 costj/),  rad 
c 

incremental  change  in  differential  collective 
pitch  angle,  rad 

incremental  change  in  collective  pitch  angle, 
rad 
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AV«J| ,. 


response  error:  measure  of  change  in  blade 

response  from  one  iteration  to  the  next* 
nondimensional 

generalized  coordinate  for  the  jth  mode  at  ith 
azimuthal  station  for  Kth  iteration 

generalized  coordinate  for  the  Jth  mode 


generalized  coordinate  of  the  Kth  mode 


0A'eB 


flatwise  slope,  (nSl  = 3v/3r) 

rigid-body  pitch  angle  along  the  blade,  rad 

angles  used  in  vortex  induced  velocity 
determination  (see  Figure  4) , rad 

cyclic  pitch  (6c  = -AjCosij,  - B^inilO  , rad 

blade  root  pitch  angle,  ( e o + e ),  rad 

o 

collective  pitch,  rad 


rigid-body  pitch  angle  plus  torsional  deflection 
at  a blade  section,  rad 

advance  ratio,  Vf/QR 

distance  from  the  raidchord,  measured  parallel  to 
the  rotor  plane,  positive  toward  trailing  edge,  ft 

2 4 

air  density,  lb  sec  /ft 

wake  influence  coefficient:  downwash  at  the 

blade  due  to  a vortex  element  in  the  wake 

average  aerodynamic  damping  coefficient  for  the 
Kth  mode 

damping  coefficient  which  couples  the  torsional 
components  of  modes  K and  J 

angle  between  chord  and  rotor  plane,  rad 

built-in  twist,  positive  for  leading  edge  up,  rad 

torsional  deflection,  positive  for  leading  edge 
up  , rad 
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I 


V 

+ 


u,. 


edgewise  slope  (V  = 3w/3r) 

azimuth  angle  measured  from  the  downstream 
position , positive  in  the  direction  of  rotation 
of  the  particular  rotor,  deg 

rotor  rotational  speed,  rad/sec 

natural  frequency  of  the  Kth  mode,  rad/sec 

indicates  time  derivative,  i.e.,  v = dv/dt 
indicates  second  time  derivative,  i.e.,  v = dv/dt 


% 


I 


Subscripts 


i 

J 

j 

K 

k 

i 

m 

n 

N 

N, 


y 

z 


denotes  radial  location 
normal  mode  number 
denotes  azimuthal  location 
normal  mode  number 
denotes  azimuthal  location 
denotes  radial  location 
denotes  radial  location 
denotes  azimuthal  location 
total  number  of  wake  elements 
denotes  number  of  azimuthal  steps 

denotes  quantities  parallel  to  the  x-axis  or 
normal  to  the  rotor  shaft 

denotes  quantities  parallel  to  the  y-axis 

denotes  quantities  parallel  to  the  z-axis 
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